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CONSERVATION LAWS FEYNMAN’S 
MODIFIED ELECTRODYNAMICS! 


Abstract 


treatment the self-energy problem, the divergence elimi- 
nated introducing convergence factor into the integral over the virtual 
photon momentum space. Feynman has remarked that his choice con- 
vergence factor inconsistent with the conservation energy for the radiation 
field atom. This problem examined more general way. The 
modification the Maxwell equations caused the convergence factor de- 
duced. The modified field equations belong the generalized electrodynamics 
described Podolsky. The modified energy-momentum tensor shown 
satisfy the conservation law for the field with source. 


Feynman’s (1, treatment the self-energy problem, the divergence 
which occurs for high momentum exchanges between electron and virtual 
photon eliminated introducing relativistic cutoff. This cutoff 
equivalent expressing the self-energy the difference between the effects 
photons zero rest mass and those rest mass divergencies then 
cancel and the self-energy the free electron can expressed change 
mass independent the momentum. 


the expression for the self-energy bound electron, the largest term, 
which depends logarithmically can eliminated recalculating the 
term energies with the renormalized mass. the remaining terms, the real 
part gives the term shift and the imaginary part gives the radiation resistance. 


Feynman has remarked, however, that this modification quantum electro- 
dynamics incomplete because the radiation resistance depends slightly 
the cutoff whereas the radiation field great distances independent the 
cutoff. Thus the energy lost atom would not check with the energy 
absorbed distant absorbers. 


This remark consists two separate points. Firstly, the source strength 
the radiating atom depends the cutoff. Secondly, the radiation field 
given source cannot conserved because depends the cutoff near the 
source, but great distances independent the cutoff. only the second 


-point that discussed this paper. 


received June 21, 1951. 
Contribution from the Department Physics, University British Columbia, Vane 
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For this point need investigate the conservation law for the field 
given arbitrary source. Feynman’s notation will used throughout. 
consider particular vector potential solution the Maxwell equations, 


which may written 


where the four-vector current source which conserved and 
the singular solution used Feynman express the interaction 
between two charges, which consists incoming waves from the past and 
outgoing waves the future. Thus A,(2) the sum the field emitted 
J,(1) the past and the field the future. 


The modification the classical electrodynamics results from the replace- 
ment the singular some nonsingular function, find the corre- 
sponding effect the field equations, simplest work momentum 
space. The Fourier transform A,, defined 


form the function and the transform the additional 47? 
comes from the integration over 


Feynman chooses for the transform for the self-energy problem the 


Multiplication the left and shows that satisfies the 
equations 


The second equation follows from the current conservation, 
These are the transforms the field equations 
(2b) 


which are satisfied given but with for 


take for the energy-momentum tensor the field the Maxwell tension 
tensor, with the modified field strengths substituted for the Maxwell field 
strengths, (the primes are dropped for the rest the paper), 
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4aTy = For + Fas Fas Sy, (3) 
then the four-divergence this tensor 


shown the appendix. 


The four-divergence the energy-momentum charge free space does 
not vanish. Thus, for source surrounded distant absorbers, energy- 
momentum would lost the space between. The term which spoils the 
conservation law vanishes allowed become infinite after mass 
renormalization. This result verifies the second point Feynman’s remark. 


However Equations (2), the field equations, belong the generalized 
electrodynamics described Podolsky and others (3). The appropriate 
energy-momentum tensor for this field has been derived Podolsky from 
the Lagrangian which gives the field equations (Equations (2)). Translated 
into Feynman’s notation, and with for this 


2 


The field then satisfies the energy-momentum conservation law 


“uy 
shown the appendix. 

The energy-momentum the field may therefore redefined way 
consistent with the modified electrodynamics that the conservation law 
satisfied. The expression for the force density, formally the same 
for the Maxwell field, but the field strengths are defined the non-Maxwell 
Equations (2). 
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Appendix 
PROOF CONSERVATION LAW 


With the field strengths defined 


where 


have the field equations 


From (1) and follow the useful relations 


Differentiating the tensor expression gives 

The remaining terms can rewritten, with the help (1), 


The first two brackets vanish (4) and the third (2) and 
conservation law then follows from (2). 


For the Maxwell tension tensor the terms are absent the four- 
divergence simply Fya Equation (4) the text then follows 
from (2) the appendix. 
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NUCLEAR MAGNETIC MOMENT SCANDIUM MASS 


Abstract 


the method nuclear magnetic resonance, the magnetic moment 
(without diamagnetic correction) found 4.74916 0.00012. The correc- 
tion 0.00717 with unknown and possibly large error. The equipment 
designed search for magnetic resonance varying the field the magnet 
described, with special emphasis the magnet current regulator, 


Introduction 


Nuclear magnetic resonance absorption (3,4) has been widely used recent 
years determine nuclear moments. The principles and general practice are 
well described the literature and will not repeated here. One the 
greatest difficulties arises from the fact that the point resonance not 
usually very well known advance, and therefore search must conducted. 
Most earlier experimenters (10,11,19) have chosen fix the magnetic field and 
vary the frequency. This inevitably results loss flexibility the radio- 
frequency apparatus. the present work magnet current regulator was 
designed and constructed permit slow, regular variation magnetic 
field over range This allowed the use radio-frequency bridge 
which would have been almost impossible tune automatically. practice 
this regulator maintained the field within one part 140,000; such system 
might well find use other fields and therefore described detail below. 


Apparatus 
Magnet 


available magnet with 334 in. poles and adjustable gap had previously 
been used work this type (1). With in. gap, required amp. 
for 8000 gauss. This power was supplied motor-generator set and 
was regulated controlling the generator field current. its best, this 
system never regulated better than one part 10,000, and was usually 
worse. This regulation was not considered good enough for automatic search- 
ing, though was entirely adequate for observing resonances the oscillo- 
scope. 


order that the magnet current might carried vacuum tube, the 
coils were rewound with number Formex wire; 500 ma. 900 then gave 
8300 gauss. practice, fields above 7000 gauss were not used, since this 


region saturation the iron the poles made the field much less uniform. 


Manuscript received April 1951. 
Contribution from the Radiation Laboratory, McGill University, Montreal, Que. 
Fellowship under the National Research Council Canada, 
address: Department Physics, University Saskatchewan, Saskatoon, Sask. 
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Exact parallelism the pole faces was found very important for uni- 
formity. They were carefully ground flat and held parallel three-legged 
spacer. Ring shims were used extend the region the uniform field. 
trial and error, the best thickness was found 1.5 times greater than those 
calculated from the article Rose (13). 


very convenient device was used check the uniformity the field. This 
was double probe consisting two coils series—placed in. apart—and 
each containing water (plus cupric chloride) sample about mm. diameter. 
This probe could rotated about axis passing through one the water 
samples. regenerative detector (1, was connected it, and this gave 
separate peak the oscilloscope for each sample, provided they were located 
different fields. Differences 0.2 gauss could detected this way. 


Magnet Current Regulator 


block diagram the regulator shown Fig. Part the voltage 
across the shunt compared with the voltage Weston standard cell; the 
difference chopped into cycle a-c. Brown converter vibrator. This 
a-c. amplified and converted back d-c. phase sensitive detector which 
preserves the original polarity. After further amplification, this voltage 
controls the magnet current. 


MAGNET 
POWER 
SUPPLY 


MAGNET 


CROSSOVER 
NETWORKS 


GAIN 5000 


ORIVER 
CONTROL TUBE 
GAIN 0-2 


GAIN= 350,000 
TIVE DETECTOR 


Fic. Block diagram the magnet current regulator. 


This system appears the only one sufficiently stable for this application 
(5). has, however, two defects: first, the d-c. output the phase sensitive 
detector contains large 120 cycle ripple component; and second, the high- 
frequency (i.e. transient) response poor. The more the ripple filtered out, 
the worse becomes the high-frequency response. Both defects are remedied 
a-c. channel parallel with the above. its regulating action, 
suppresses the ripple factor 5000, and the same time supplies the 
high-frequency response. 

The loop gain 100,000 for d-c. and 10,000 for a-c. Thus, output fluctua- 
tions will reduced these factors. For 10% line voltage change, the magnet 
current should change one part million; rapid changes should 
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have the same effect. ripple from the phase-sensitive detector held 
the same amount. the voltage (on the shunt, Fig. considered the 
output, the feedback factor for a-c. and d-c. 0.02 (when the potentio- 
meter center-scale). This leads the gains shown the diagram for the 
various amplifiers, since the Brown sensitive detector com- 
bination has gain actual practice, the gains were adjusted give 
best operation and were not measured; hence the values shown are the design 
values only. 


The schematic circuit given Fig. with the blocks Fig. identified. 
Because the wide range currents handled single control (i.e. 
without change supply voltage), the control tube must have large plate 
dissipation rating; hence tubes such are not suitable. 304TH 
very satisfactory, but the higher currents the grid must driven about 
positive and then draws about ma. This supplied the 6L6 cathode 
follower. The 180 ohm cathode resistor allows the grid negative while 
staying positive with respect ground. The potentiometer the grid cir- 
cuit used cut off the magnet current slowly, order avoid dangerous 
switching transients. When turned all the way up, cam the shaft 
closes microswitch, that the sliding contact need not carry the grid current. 


The shunt wound with insulated manganin wire water cooled brass 
tube. tap ohms feeds the d-c. comparing circuit, while the full 100 ohms 
feeds the a-c. channel. The network parallel with the magnet makes the load 
constant resistance 1790 ohms all frequencies, and produces smooth 
and nearly flat response 20,000 cycles, the limit the measurements. 
Without the network, the gain varies with frequency extremely compli- 
cated way, evidently result self-resonances the magnet. Since rapid 
changes current are not expected arise the magnet itself, the regulation 
the magnet current should not affected. 


The potentiometer the comparing circuit turn Helipot, driven 
clock motor through gear train and friction clutch which allows manual 
adjustments. The scanning rate may changed changing gears and 
switching additional resistance series. These wire-wound precision resistors 
also reduce the effect the steps the Helipot, which are quite perceptible 
when used alone. 


The cycle amplifier has passband from 700 cycles; high frequencies 
are cut off keep stray voltages from affecting the detector. This simply 
grid-controlled full-wave rectifier; the plate voltage supplied small 
receiver power transformer. The small phase shift that was required pro- 
duced the network the primary circuit. The variable resistor the 
cathode circuit serves manual current control; for best operation should 
set (with the regulator not operating) give current near that used. 


The a-c. amplifier has pass band from 0.1 200 cycles; the high-frequency 
cutoff goes average db. per octave about one megacycle. This has 
important bearing the stability the regulator and will discussed 


below. 
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order avoid trouble with hum, which would particularly bad the 
cycle amplifier, the five 12SH7’s are heated with d-c. from the high-voltage 
supply. 12SH7’s were used because they happened plentiful; but, any 
case, the high-frequency cutoff the a-c. channel requires tubes high trans- 


Because the unbypassed cathode resistor, the mixer stage adds the volt- 
ages applied the grids. The crossover networks, shown enclosed dotted 
lines, shape the cutoffs the two channels give smooth transition over- 
all gain from 100,000 10,000. This important since improper choice 
filter constants could easily produce enough phase shift, opposite sign the 
two channels, give mutual cancellation gain and drop the over-all 
characteristic below the gain either channel alone. crossover fre- 
quency about 0.5 cycle per sec.; thus the low frequency filter removes good 
deal the 120 cycle ripple from the phase sensitive detector, and the a-c. 
channel removes the rest. 


order that the regulator should not oscillate, the loop gain should less 
than any frequency where the phase shift 180°. high-frequency cutoff 
db. per octave will produce such phase shift; any smaller slope safe. 
already noted, the actual slope db. per octave one megacycle; the 
loop gain down 100 kc. 


The d-c. power requirements are follows: for the magnet, 1500 v., 500 
for the regulator, 350 v., 250 ma. (6L6 and heaters), 300 (regulated) 
ma., 150 (VR tube) ma., and 150 (VR tube) The 
bias voltages and 3.9 are derived from the last; the 3.9 tap had 
decoupled 200 uf. condenser ground. 


Observation the proton resonance the dscilloscope showed that the 
magnetic field variations were about 1/20 gauss peak-to-peak, out 7000 gauss 
total. Long-term stability was probably not quite good; but any case, 
was not measured since was importance for this application. Scanning 
was smooth and even, and readings could reproduced accurately the 
dial could read. 


This system would well adapted control the proton resonance. 
was felt, however, that scanning would then good deal more difficult. 
likely that the radio-frequency spectrograph Pound (10) would work equally 
well. 


Other Apparatus 


twin-T bridge (2) was used; was found very much easier adjust 
than the type described Bloembergen, Purcell, and Pound (4). The radio- 


frequency source was General Radio 1001-A signal generator. This was fairly 


satisfactory, but did not have good frequency stability could desired. 
low-noise preamplifier (18) was followed Hallicrafters SX-42 receiver, 
lock-in amplifier, and recording meter the usual arrangement. The magnetic 
field was modulated cycles small coils around the magnet poles. 
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The regenerative detector already mentioned was used, with different probe, 
give the proton resonance for measurement magnetic moments. Since the 
field was slightly different the unknown sample, the proton resonance was 
observed simultaneously with the bridge. The regenerative detector could 
made lock with the signal generator, and the difference magnetic field 
measured the oscilloscope. never exceeded 0.5 gauss. 


Method Search 


Since magnetic moments are known advance only +10%, and since 
one must not sweep over resonance less than several times the time con- 
stant the indicating circuit, long searches are necessary. Experiments the 
proton resonance showed that could considerably broadened, with 
visible loss amplitude the recorder, using modulation amplitudes 
than the line width; gauss was about the best value. scanning rate 
gauss per four time constants then gave trace which was distorted but 
still easily recognized. This improvement least times over the 
conditions that give undistorted record. 


The big advantage bridge over the circuits that scan frequency (10, 
11) that can observe the dispersion component resonance. Thus the 
radio-frequency voltage applied need only chosen large get the largest 
possible signal. For absorption, the sensitivity best only for certain voltage 
which not known advance. the other hand, the bridge balance not 
stable enough allow left unattended. 


Experiments Sc* 


Scandium has one isotope mass 45. Spectroscopic measurements (7, 14) 
give spin 7/2 and magnetic moment 4.8 nuclear magnetons. Theory (4) 
thus indicates that the signal should strong. solution scandium nitrate 
was prepared dissolving scandium oxide (spectroscopically pure) nitric 
acid. The resonance was easily found near the expected place. Comparison 
with the proton resonance described below gave the frequency ratios shown 


TABLE 
OBSERVATIONS THE SCANDIUM RESONANCE 


Run No. Scandium frequency, Proton frequency, Ratio Sc/P 
Mc. Mc. 
6.90100 28.40740 0.242930 
6.90000 28.40356 0.242927 
6.89930 28.40128 0.242922 
6.89930 28.40141 0.242943 
6.89929 28.40124 0.242945 
6.89930 28.40127 0.242944 
6.89932 28.40156 0.242951 
6.89933 28.40090 0.242948 


Mean ratio 0.242939 0.000003 
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Table Measurements were taken rapid succession, each one consuming 
about min. There was evidence drift the magnetic field. 


The Helipot dial setting for the center the scandium resonance was noted 
the field was slowly scanned past; this was done three times make sure 
that everything was stable. The dial was set this point and the proton 
resonance centered gauss oscilloscope sweep, and the two frequencies 
measured. 

the g-factor the proton taken (16) 5.58536 0.00012, the magnetic 
moment 4.74916 0.00012 nuclear magnetons. The diamagnetic 
correction (9, 17) 0.00717, with probable error difficult assign. 
taken 10%, the final result 4.756 0.0007 nuclear magnetons. The 
spectroscopic value, 4.8, thus only low. 


Discussion 


view the recently observed dependence resonance frequency 
chemical compound, this result not accurate indicated the probable 
error observation. not possible present say how large such 
error might be. 

Other workers have independently found this resonance and measured the 
magnetic moment. Proctor and (12) find frequency ratio 0.24289 
0.00004, with scandium nitrate, and Sheriff and Williams (15) 0.24296 
0.00005 with the chloride. Both agree within their errors with the more 
accurate measurements reported here. 

The author indebted many members the Radiation Laboratory for 
stimulating discussions, the machine shop staff for their careful work, and 
Dr. Thompson for help with the sample. wishes thank especially 
Dr. Anderson, who built the original apparatus, Mr. Bloom, who 
helped with the rebuilding, and Dr. Foster, who supervised the work and 
contributed valuable suggestions and encouragement. 
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FUNDAMENTAL THEORY THE MAGNETISM 
MASSIVE ROTATING 


GEORGE LUCHAK? 


Abstract 


phenomenological theory, based relativistically covariant generalization 
equations include gravitational fields, developed account for 
the magnetic fields massive rotating bodies. The equations yield the Wilson- 
Blackett expression for the magnetic moment the earth and stars but give 
magnetic field for mass-bodies moving without rotation their own gravitation- 
fields. They indicate that the magnetic field due the motion the earth 
its orbit negligibly small compared the field due its rotational motion, 
and they provide possible explanation for the variable magnetic fields 
light-variable stars. 


Blackett (2) has recently drawn attention the approximate validity for 
the earth, sun, and five stars, the empirical formula 


between the magnetic moment and the angular where 
the gravitational constant, the velocity light, and constant 
the order unity characteristic the body. and range relative 
value from and the magnetic fields, which the values are 
based, range relative value from 1800 From detailed discussion, 
Blackett concluded that Equation (1) was the manifestation new property 
matter not yet included present day physical theory, and, since con- 
tained the gravitational constant, might well the result interaction 
between gravitational and electromagnetic fields. 


simple hypothesis, from which Equation (1) can derived, was given 
Wilson (14) who postulated that mass element moving with velocity 
produced magnetic field distance given 


analogy with the magnetic field moving charge. Wilson’s equation 
can also written the equivalent form 


where the mass density, the velocity, and fictitious current due 
the motion the mass. Equation (3) and, through it, Equation (2), are 
not satisfactory physical laws because they are not relativistically covariant, 
and predict magnetic field due translational motion. Swann (10) tried 
remedy the objections Wilson’s theory but his subsequent equations 
proved physically untenable The most recent attempt improve 
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the Wilson expression was made Papapetrou (7) who postulates that 
every mass element there corresponds unique point which 
gravitational center. Then, the four dimensional velocity vector 
the mass element, the corresponding vector for the gravitational center 
and quantity satisfying the equation 


(4) 


where the metric tensor flat space and the four dimensional 
current vector replace the three dimensional current vector Equation 
can shown (6) that Equations (4) (6) are not relativistically 
covariant and, furthermore, not determine uniquely for any given physical 
situation. 


Since there appears evidence the quantum mechanics materi- 
als ordinary temperatures and pressures for interaction which might 
lead Equation (1), probable that such effect manifests itself only 
under the extreme physical conditions existing the earth’s core and the 
sun and stars. Ramsey (8) has recently shown that, the boundary the 
core, the effect high pressure sufficient change radically the properties 
matter from molecular metallic state. The presence such tran- 
sition might well lead quantum effects analogous those that occur 
low temperatures. Since, however, suitable quantum mechanical theory 
superconductivity not yet available, seems unlikely that quantum 
theory explanation for Equation (1) found easily. possible, 
the other hand, proceed from the other end and produce phenomeno- 
logical theory, was done London and London (5) describe supercon- 
ductivity, suitable generalization equations. Into these 
equations will necessary incorporate gravitational field relativis- 
tically covariant fashion and, addition, interaction the proper sort 
between the electromagnetic and gravitational vectors. One consequence 
such interaction, according Tzu (12), the necessity introducing 
least one new parameter addition and Referring again the 
London theory, interesting note that the extension Maxwell’s equa- 
tions superconducting media accompanied the introduction 
parameter whose value depends the substance and its physical state. 
here suggested that the new parameter which appear the inter- 


_action term theory which might lead Equation (1) not universal 


constant for all matter all states, but varies with physical conditions and 
only appreciably different from zero under conditions similar those that 
obtain the sun, stars, and the core the earth. This contrary the 
theories Wilson, Swann, and Papapetrou. Recent work Runcorn (9) 
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demonstrating that the outer crust the earth does not contribute the 
main geomagnetic field corroborating evidence. 


Formulation the Equations 


Maxwell’s equations for macroscopic media, expressed tensor form, are 


+ + Feed 0, (7) 
these are added 
Ke ad (9) 


—B, 0 B, —ik, 
—B, 0 —tk, 


ab — z z y 


Far (10) 


where imaginary time ict used, where the metric tensor 


flat space and the force density, the charge density, the component 
the current density, the component the electric field, the compo- 
nent the magnetic field, the component the magnetic induction, 
and the component the electric displacement vector when Equations 
(7), (8), and (9) are written out Cartesian coordinates, and gaussian units 
are used. formal distinction between covariant and contravariant tensors 
unnecessary for flat space but retained order that possible extension 
the equations general covariance will obvious. 


Following Corben (3), Equations (7), (8), and (9) are generalized include 
gravitational fields follows. They are extended five dimensions while 
are retained. the extension, and are kept skew symmetric. Thus, 
obvious that and are four-vectors and invariant when the 
fifth dimension transforms into itself and the four physical dimensions trans- 
form Lorentz space. Maxwell’s equations five-space are 
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When these equations are rewritten that five-dependence the tensors 
coordinates shown explicitly, get 


Equation (9) written out five dimensions the familiar vector form, 

and W=J-E+ (22) 


where three-vector whose components are and Thus 
would appear the gravitational force per unit volume and 
the rate which gravitational forces work unit volume. This would 


suggest that four-vector 


introduced, where the rest-mass density e.s.u., and thus invariant, 
and the four dimensional gravitational field vector. Equation (20) 
reduces four dimensional version the Poisson equation 


(24) 
Writing Equations (23) and (24) become 
where invariant. worth noting that the relation 
(27) 


analogous the ‘pair’ properties and for isotropic media 
where and are scalar functions position which characterize the electro- 
magnetic properties matter the large. The inference drawn that 
still another scalar the same category. For simplicity will consider, 
all that follows, only homogeneous media for which constant. The 
appropriate generalization for nonisotropic, nonhomogeneous media will 
obvious. Thus Equations (17) (20) become 
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generalization Maxwell’s equations with gravitational forces but with 
interaction between the electromagnetic and gravitational fields. Such 
interpretation inadequate for the purpose hand. The simplest remaining 
interpretation take either identically zero and choose the 
other suitable manner. Suppose Then Equations (28) and 
(30) become 


and 


These are the ordinary Maxwell equations. Thus, there are impressed 
charges currents, the electromagnetic field will zero. Such theory 
could not account for the expression other alterna- 
there interaction, must depend §*. The dimensions must 
and, according the expression, must vanish 
identically for nonrotating sphere the time-independent state. Thus, the 
interaction term must depend the kinematic conditions. The simplest 
possible kinematic vector the four-velocity and assume that appears 
the interaction term the simplest possible way. arrive inter- 
action which also four-vector, must multiplied either scalar 
which involves tensor rank two formed from The simplest 
scalar that can constructed This leads four dimen- 
sional Wilson—Blackett expression which rejected unsatisfactory (7). The 
simplest tensor rank two that can formed from The expression 
fulfills all the required conditions and Equations (28) (31) can now 
written their final form 


where all indices run from 


Equations (32) (35) are obviously internally consistent, since (33) and 
(35) can solved uniquely, shown Corben (3), and (32) and (34) 
are merely Maxwell’s equations for moving media with additional four 
dimensional current The four dimensional divergence the latter 
not necessarily zero. consequence that, the above description, 
charges may induced the motion matter through gravitational 
fields. equations for moving media, without gravitational inter- 
action, exhibit the same phenomenon. 


Equations (32) (35), written the familiar vector form, become 


— 
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(37) 

(41) 


where (V;, V,, the three dimensional velocity vector, terms the 
order V/c|? greater have been neglected, and 


(43) 


the well known three dimensional current density for moving media, where 
the electrical conductivity. 


Since Equations (36) (42) are covariant for special relativity, im- 
portant designate which are inertial systems. These are chosen, 
our case, from among those that have inherent angular velocity. Since 
allowable speak velocities rotation absolute manner, without 
the necessity specifying any framework with respect which the angular 
displacement takes place, the distinction valid (13). 


II. Magnetic Field Sphere Rotating its own Gravitational Field 


Consider sphere for which whose angular velocity constant, 
and for which the rest-mass density constant within the sphere and zero 
outside. Assume that conditions are independent time. Equations (36) 
(42), for observer stationary the center the sphere, then become 


. 
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and curl (50) 


where wy, wx, 0), the rest-mass density c.g.s units 


Equations (49) and (50) can solved give which, when 


substituted Equation (46), 


and div (54) 


From Equation (53) can seen, taking the divergence both sides, 
that 


also evident that 


where [B,, the vector spherical coordinates and there symmetry 

Equations (55) and (56) have the solution 


(57) 


Thus Equations (51) (54) reduce 


and div (61) 


| 
| 
| 
| 
| 
| 
| | 
| 
| 
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Equations (58) (61) lead the Wilson—Blackett expression (1). 


Assuming that the core the earth and the crust and 
the atmosphere, turns out that the core. For the sun, assuming 
polar magnetic field gauss, 1.6, and for the average for five stars 
Blackett (2)). for the earth and sun, seems 
approximately numerically equal the density for material high tempera- 
tures and pressures. 


III. The Magnetic Field Due the Motion the Earth 
the Gravitational Field 


The earth, its motion about the sun, will generate gravitational waves (3). 
first approximation, however, neglecting terms (V/c) where the 
velocity mass element the earth and the velocity light, the gravi- 
tational field obtained taking the static field the earth and transporting 
along with the earth its motion. The gravitational field acting any 
mass element now composed two parts; the gravitational field due the 
sun, and the gravitational field due the remainder the earth. the mass 
the sun its average density and its radius then the gravitational 


MGR 
units, 


The gravitational field mass element the core due the earth itself 


where the distance from the center the earth the mass element and 
the density the core (assumed uniform). now necessary 


calculate 


earth rotating with angular velocity the velocity any point the core 
(66) 


where the velocity the mass center and and are measured with 
respect coordinate system whose origin moves with the center the 
earth and whose coordinate axes are always parallel the axes the non- 


j 
| 
| 
| 
| 
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rotating observer the center the sun. elementary calculation then 
shows that (see appendix) 


The first term Equation (67) leads the ‘rotational’ magnetic field. Taking 
1.2 for half the rotational velocity the core’s equator, then 

Thus the nonrotational components the earth’s field due interaction 
with the sun’s gravitational field are negligibly small compared the rota- 
tional component. 


2.5 


IV. Other Effects 


double star with highly elliptic orbits would exhibit periodic magnetic 
field. can seen from Equation (67) that there would times when the 
terms involving (A/R)* would small and other times when they would 
make significant contribution the magnetic field. the rotation one 
the components the double star were retrograde with respect its 
direction revolution about the other component, might possible, 
addition, observe reversal magnetic field. 


star moves without rotation its own gravitational field, the change 
vector and scalar functions following the motion zero and thus, Equa- 


tions (38) and (39), 


Such star would exhibit electromagnetic field. 
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Appendix 
DERIVATION EQuaTION (67) 
Equations (63) and (64) give 


the coordinates the mass-center the earth are given 


[Xi(t), Y(t), Zi(¢)I, 


Thus Equation (67) readily follows. 


1 
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THE HALF-LIFE IRIDIUM 


KASTNER 


Abstract 


Lauritsen electroscope was used compare the intensity source 
with radium standard. The result least squares fitting the decay curve 
was half-life 74.37 0.07 days. 


Introduction 


radiography (6), wide deviations were noticed the various values the 
half-life quoted the literature (Table I). 


TABLE 
VARIOUS VALUES THE HALF-LIFE 


Value, days Author Date 
Fermi and Amaldi (2) 1936 
Jaeckel (5) 1938 
70+1 Friedlander al. (3) 1944 
Goodman and Pool (4) 1947 


was felt that more precise knowledge the half-life would worth- 
while view the potential industrial and medical importance this isotope. 
Since the commencement this work two more values not quite agreement 
were reported: 74.5 0.7 days Sinclair and Holloway (7) and days 
Cork al. (1). 


TABLE 
ANALYSIS 
Relative Logarithm Time, 
(relative days Residuals 
activity) 
1.834 0.2634 0.0008 
1.412 0.1498 28.01 0.0006 
1.125 0.0512 51.13 0.0045 
1.025 0.0107 62.19 0.0001 
0.942 0.0262 72.01 0.0027 
0.812 0.0904 86.02 0.0048 
0.718 0.1442 100.16 0.0013 
0.636 0.1965 114.18 0.0031 
0.558 0.2532 128.14 0.0029 
0.486 0.3131 142.14 0.0003 
0.431 0.3658 156.00 0.0030 
0.381 0.4190 170.18 0.0069 
0.327 0.4849 184.06 0.0025 


Manuscript received May 22, 1951. 
Contribution from the Division Physics, National Research Laboratories, Ottawa. 
Issued N.R.C. No. 2540. 
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Lauritsen electroscope was used compare the intensity the 
source used the radiographic experiments with that 100 mgm. radium 
standard. careful comparison was made about two-week intervals for 
observation period extending over more than six months (two and one-half 
half-lives). 


The linearity the decay curve, and thus the radioactive purity the 
source, were established the method least squares (8) from the data 
obtained. These data are shown Table II. The least squares fitting sup- 
plied the decay equation, log (relative activity) 0.2626 0.004048 days, 
from which was calculated the half-life and its probable error. The result was 
half-life 74.37 0.07 days, which agrees with and rather more 
precise than the value 74.5 0.7 days published during the course the 
present work (7). 
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SOME QUANTITATIVE MEASUREMENTS 
THREE-CENTIMETER RADAR ECHOES FROM FALLING 


Abstract 


the winter 1948-49 series measurements was made Ottawa 
attempt verify, experiment, the theoretical computations Ryde 
for the intensity backscattering cm. radio waves snowflakes. The 
radiation backscattered from falling snow was measured, together with the actual 
mass snow per cubic meter atmosphere. Whenever possible the size dis- 
tribution the snowflakes was determined. Owing the small number snow- 
storms studied, firm conclusions cannot bedrawn. However, the results obtained 
are agreement with those predicted Ryde, and provide useful practical 
information. 


Introduction 


Rainstorms cause echoes appear radar sets, operating centimeter 
wave lengths, owing radiation being backscattered from the raindrops. 
Similarly, radiation backscattered from snowflakes will cause radar echoes. 


From theoretical considerations Ryde (5) has given quantitative measures 
the above mentioned phenomena terms the radio frequency power 
backscattered from the respective particles. theoretical expression for the 
case rain, and the modifications this expression for the case snow are 
discussed Section 


Section description given experiment which was performed 
measure the power radar signals wave length 3.3 cm. returned from 
falling snow. The comparison the experimental results with those obtained 
theoretically given Section 


Theoretical Considerations 


Ryde’s theoretical expression for the power the backscattered 
radiation from rain has been modified Marshall, Langille, and Palmer (4). 
Assuming that the cross-sectional area the radar beam uniformly filled 
with falling rain and that the size the raindrops small compared with the 
wave length, the appropriate expression 

where 


backscattered power received the radar set, 
peak output power the radar, 


received June 13, 1951. 
Contribution from Defence Research Telecommunications Establishment (Radio Physics 
Laboratory), Ottawa, Ont. 
Defence Research Telecommunications Establishment (Radio Physics Laboratory). 
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effective area the antenna, 

pulse length, 

refractive index water, 

range the rain from the radar set, 

wave length, 

sum the sixth powers the diameters the raindrops 
unit volume atmosphere. 


According Ryde the expression for rain should also hold the case 
snow provided: 

(a) considered rather than where the sum the 
squares the masses unit volume atmosphere. Thus replaced 
where the density ice. 

the refractive index for ice substituted for 

(c) dimensionless factor introduced account for the nonspherical 
shape the snowflakes. According Ryde equal unity for spheres, 
and approximately two for plates and needles. (In this paper has been 
assumed equal unity). 


After collecting constants and substituting for the value 1.75 (3), and for 
the value 0.916 gm./cc., the expression the case power reflected from 


snow becomes 
Mr? 


For the calculation the following two methods may used: 
(a) assuming the snowflakes equal mass, then 


where 
total number snowflakes per unit volume atmosphere, 


total mass snow per unit volume atmosphere. 

(b) assuming that the masses the snowflakes are nonuniform and 
obtaining the snowflake mass distribution from collected samples snow- 
flakes. 

Both methods require measurements the rate snowfall, and know- 
ledge the velocity fall the snowflakes. 


Experimental 


The measurements described this report were made Ottawa, Canada, 
the winter 1948-49. modified naval radar equipment (type SU) was 
used. 


Power measurements were made comparing, A-scope presentation, 
the signal received from falling snow signal from cm. signal generator, 
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fed into the radar means wave guide directional coupler. The attenu- 
ator the signal generator was adjusted two-minute intervals match the 
amplitude the reference signal with that the signal reflected from the 
snow. The range was limited 1000 yd. order study the physical 
properties the snow. 


The physical properties the snow producing the echoes were measured 
from mobile shelter mounted truck. order clear Jocal obstructions 
was necessary elevate the radar beam 15° above the horizontal. this 
way the radar echoes were received from snow located approximately 800 ft. 
above the ground surface range 1000 yd. The mobile shelter was 
placed position such that would intercept the snow that produced the 
snow signal the radar. the calculations allowance was made for the 
time necessary for the snow reach the ground. 


The rate snowfall was determined weighing the snow that accumulated 
rubber ground sheet. The ground sheet measured in. and was 
placed flat the ground well removed from obstructions. Every five minutes 
the sheet was lifted and the snow was placed aluminum bowl. The bowl 
and the snow were weighed beam balance inside the mobile shelter. This 
measurement the rate snowfall appeared quite satisfactory for low 
wind speeds. ‘At wind speeds greater than miles per hour doubtful 
whether this method gave reliable measure the rate snowfall. Attempts 
were made develop type snow measuring device which would draw 
air and snowflakes and separate out the snow. Four models were tried and, 
although the last model was quite promising, insufficient measurements were 
made trust the results. 


The size distribution the snowflakes was measured means the dyed 
filter paper technique used for determining the size raindrops 


Samples snowflakes were caught filter papers which had been lightly 
dusted with gentian violet dye. The papers were then held over electric 
heater for moment allow the snowflakes melt, producing stains the areas 
which were proportional the masses the original snowflakes. 


For each the measurements P,, the value was determined from 
measurement the rate snowfall, over the corresponding time interval 
and from the data obtained filter paper exposed during portion this 
time interval. The snowflake stains the filter paper were divided into 


class intervals having the means the squares the 
masses the respective intervals. Thus 


are the numbers snowflakes each respective class. The subscript 
refers the filter paper measurement. 

Assuming that the filter paper represents the true snowflake distribution 
over the time interval the corresponding measurement, then 
where numerical constant which is, fact, the ratio that the 
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sample bears unit volume atmosphere. This numerical constant could 
obtained considering the size the filter paper and the length time 
exposure. However, was considered that the measurement the latter 
would constitute possible large source error the results. Hence, was 
decided obtain the value the numerical constant assuming that the 
velocity fall, the snowflakes was constant(2), and using the corres- 
ponding measurement the rate snowfall, obtain =m, the total mass 
snow per cubic meter atmosphere. 


=m = 
where given gm./cu.m. (atmosphere) 


and the equation the case snow becomes 


Correcting for units and substituting the values for the radar set parameters, 
taking m./sec. the equation becomes 


2 


mgm. 


Results 


During the winter 1948-49 measurements were made snowstorms 
but sufficient data for analysis was obtained only four storms. these 
four storms, trustworthy size distribution data were obtained only one, 


the 22nd February. 
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VOL. 


JANUARY 1949 


Decibels one milliwatt 


Fic. The data for Jan. show the power the snow echo received, Pr, plotted against the 
rate snowfall, The solid line the least squares line through the experimental points. 
all the snowflakes were equal mass, the predicted shown the brokenlines. Broken 
lines for three possible flake masses are shown. this day the approximate temperature 800 
ft. altitude was 14°C. 


JANUARY 1949 


(as water) 


Obs 80 75 70 65 60 
Decibels below one milliwatt 
Fic. The data for Jan. are similar those number beside some experi- 


mental points indicates more than one point this position. The temperature 800 ft. altitude 
was about 10°C. 
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Fig. shows the results for Jan. 21. this day insufficient information 
regarding size distribution was obtained for the purposes comparison with 
theory and was assumed that all the snowflakes were equal mass. Mass 
sizes were chosen agree with the experimental results. The masses required 
are reasonable order magnitude but the experimental] results require 
that the mass the flakes should increase the rate snowfall increases. 


Fig. displays the results for Jan. 31. The situation was similar that 
Jan. except that was necessary assume flakes greater mass fit 
the experimental results. this day also, the experimental results require 
that the mass the flakes should increase the rate snowfall increases. 


FEBRUARY 22, 


Decibel below one milliwatt 


Fic. This shows all experimental points for Feb. 22. The points for which size distributions 
were obtained are shown Fig. The temperature 800 ft. was approximately 10°C. 
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Fig. shows all the results for Feb. The correlation between and 
poor this day and least squares line has been drawn. 


Fig. shows the results for Mar. and The slope the least squares 
line different from that Jan. and and indicates that the masses the 
snowflakes were more uniform than the other days, and according theory 
should about .05 mgm. 


Size distribution data were available for only portion the power measure- 
ments for Feb. 22. These measurements have been plotted Fig. without 
applying the correction for the size distribution data. was found that 
and gave correlation coefficient 0.67. Fig. the points Fig. 
have been corrected for the size distribution the snowflakes, and and 

four decibels between the experimental and theoretical results and possible 
that this may have been due part the method measurement which 


gave correlation coefficient 0.79. The result predicted theory 


Decibels below one 
Fic. The data for March and are similar those Jan. and 31. The number 
beside some points indicates more than one point this position. The temperature 800 ft. was 
approximately 
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pulsed reference signal was compared with noise-type snow signal. 
interesting note that Hood (1) measurements the radar echo inten- 
sities rain and cm. found the experimental results low 4.5 db. 


comparison rain and snow signals was attempted using both horizontally 
and circularly polarized radiation cm. wave Unfortunately the 
quarter-wave plate was not constructed until very late the winter season 
and only few results were obtained. The table below characteristic 
the results obtained. 


Horizontal Circular 
polarization polarization 


Snow echo db. db. 
(above noise) 
Rain echo db. db. 


FEBRUARY 22, 1949 


(as water) 


Decibels below one milliwatt 


Fic. above points have not been corrected for size distribution data. The line the least 
squares line for the experimental points. The correlation coefficient for and 0.67. 
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FEBRUARY 22, 1949 


0.05 


Decibels below one milliwatt 


Fic. This graph shows the effect applying the size distribution correction the points 
The broken line the theoretical result predicted Ryde. The correlation coefficient for 
and 0.79. 


Conclusions 


Although insufficient snowstorms were studied for firm conclusions, the 
results are agreement magnitude with those predicted Ryde. The 
size distribution snowflakes appears important calculating the radar 
echo expected from particular snowstorm. However, reasonable 
estimates may made from rates snowfall alone. would seem that the 
snow echo for any rate snowfall for any particular radar set could cal- 
culated with fair accuracy employing theoretical expression. 
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SHADOW-CASTING WITH HALIDES THE ELECTRON 
MICROSCOPE! 


BARNES AND MARGARET LAMBE? 


Abstract 


Experiments the behavior halides shadow-casting agents electron 
microscopy are described. The use calcium fluoride, placed directly the 
screen, suggested convenient method for routine qualitative 
shadowing preparations the preliminary stages electron microscopic 
study. vacuum coating equipment required. 


Introduction 


The observation that certain substances the electron microscope show 
autoshadowing not new (2). The effect most noticeable under intense 
electron bombardment from biased gun. The present investigation arose from 
routine examination cryptocrystalline aggregates, interspersed with long 
slender needles, lithium fluoride obtained from aqueous solution Dr. 
Stedman the Applied Physical Chemistry Laboratory, National Research 
Council. Heavy shadowing all samples was observed the electron micro- 
scope. was traced immediately relatively large masses the lithium 
fluoride selecting several areas the specimen screens and following the 
shadows back their origins appropriate movement the specimen stage. 


The exceptionally strong autoshadowing lithium fluoride, which 
example shown Fig. (top left), suggested that this salt might em- 
ployed shadowing agent for other objects. was realized, course, that 
the effect would only qualitative since there could control over the 
angle shadowing. Owing the simplicity adding one two particles 
salt directly specimen screen, however, compared with the usual 
vacuum coating equipment and shadowing procedure, seemed worthwhile 
investigate this possibility more detail and examine the shadowing 
possibilities other halides. 


Experimental 


Dow polystyrene latex particles (580G, Lot 3584) (1), were selected objects 
for the tests. The following halides were chosen potential shadowing agents 
largely the basis ready availability reagent, reasonably pure, grades: 
fluorides aluminum, barium, cadmium, calcium, chromium, lithium, potas- 
sium, sodium, zinc; chlorides calcium, cadmium, chromium, cobalt, copper, 
potassium, silver, sodium; bromides cadmium, copper, lead, mercury, potas- 
sium, sodium; iodides barium, cadmium, calcium, sodium. 


Manuscript received July 1951. 
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The individual salts were added dry the specimen screens the form 
fine powders very small crystals. Best results were obtained the use 
very small quantities (0.05 mgm. less) the form individual crystal 
aggregate. the latter case, shadowing was not successful unless the 
size the constituent particles was larger than about two microns. The salt 
was placed the Formvar supporting film covering the specimen screen after 
drop the latex suspension had been added and allowed dry. Upon inser- 
tion the specimen the microscope the intensity the electron beam was 
kept low possible during initial focusing. Some the salts shadowed even 
under these conditions; with others, shadowing did not occur until the beam 
intensity was increased. 


Although shadowing was observed individual screens most the salts 
examined, only very few the halides gave consistently good results. 


general, the fluorides shadowed most easily. They usually gave reasonably 
even and fine grained backgrounds except when too large particle was em- 
ployed with resultant excessive deposit the objects and supporting film. 
Sodium and calcium fluorides were consistently the best. 


was difficult obtain appreciable shadowing with the chlorides even 
under intense electron bombardment. When shadowing did occur the back- 
ground frequently was very grainy and sometimes coarsely crystalline. Very 
fine powders usually gave shadowing while those having larger particle 
size tended cause the electron beam behave erratically with consequent 
rapid variations intensity and focus. Both sodium and calcium chlorides 
appeared have direct effect the latex particles causing them become 
fluid, spread out smears over the screen, and coalesce. 


The bromides and iodides shadowed more easily than the chlorides although 
rather intense electron beams were required most cases. The quality the 
shadowed images, however, generally was very poor owing heavy, granular, 
and uneven backgrounds. 


For particular anion, best results were obtained with salts the alkali 
and alkaline earth metals. For particular cation, the fluoride was the most 
effective salt. 


few examples the shadowing effect are presented Fig. (fluorides) 
and Fig. (chlorides, bromides, iodides). The magnification factor for each 
photograph indicated line, representing length approximately one 
micron, below each picture. The reproductions are the form negative 
prints. The shadowed objects are Dow latex particles each case except Fig. 
top left, which illustrates the autoshadowing lithium fluoride, and Fig. 
second from bottom right, for which bacterium coli, four hour culture) 
was employed. 


The photographs have been selected from very large number micro- 
graphs representing typical cross section the more satisfactory obser- 


5 
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Shadow-casting with fluorides. 


PLATE 
| 
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Fic. Shadow-casting with chlorides, bromides, and iodides. 
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Thus, Fig. the occasional appearance shadowing trom two 
directions illustrated with sodium fluoride and barium fluoride. Excessive 
coating the latex particles and background when too large quantity the 
shadow-casting agent used shown with calcium fluoride. markedly flat 
field with doubled shadows (possibly artifact) appears Fig. sodium 
iodide. similar effect was observed one two photographs (not shown) 
with other iodides. other occasions, however, the same salts gave better 
results (e.g., Fig. barium iodide). 


Conclusion 


The shadow-casting operation most successful the case the fluorides 
probably because their relatively high thermal stability which ensures simple 
sublimation the salts without decomposition. That shadowing does not 
occur more than one two directions given screen indicates that only 
relatively high peaks the shadow-casting agent are involved. Scattering and 
absorption electrons the comparatively large mass the salt gives mole- 
cules the surface sufficient kinetic energy escape faster rate than they 
return. the point departure relatively high and sharp many these 
molecules, travelling straight lines, reach the surface the screen 
oblique angle, transfer kinetic energy the surface, and are trapped surface 
forces faster rate than they can escape. 


The poor quality the shadowing bromides and iodides probably due 
the lower stability these salts. fact, observations the electron micro- 
scope indicate that sublimation accompanied (probably preceded) de- 
composition that the shadowing due large measure the metal which 
some cases (e.g., mercury, lead) appeared molten the screens. 


The behavior the chlorides somewhat surprising view the relative 
ease with which they can sublimed, without decomposition, heating 
air. possible that better results might obtained with specimens other 


than the Dow latex. 


From practical point view the only salts examined that are worth con- 
sideration shadow-casting agents are lithium fluoride, sodium fluoride, and 
calcium fluoride. Calcium fluoride probably the most readily available the 
average laboratory and generally gives higher proportion satisfactory 
micrographs than does the lithium salt. For routine use has the great advan- 
tage over both lithium and sodium fluorides much lower solubility, and 
hence toxicity. Calcium fluoride has been adopted for use this laboratory 
convenient means enhancing contrast the preliminary stages 
investigation. would useful laboratories having facilities for orthodox 
shadow-casting procedures. vacuum coating unit required and virtually 
time consumed the shadow-casting process. contamination the 
microscope attributable the use these salts has been observed date. 
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IONIZING POWER COSMIC RAY PARTICLES SEA 


CHATTERJEE? 


Abstract 


Using proportional counter telescope arrangement, experiments have been 
carried out sea level explore the nature and ionizing power particles 
the soft component cosmic radiation and those produced under 1.8 cm. and 
cm. lead. The results indicate preponderance relativistic electrons 
the soft component and under cm. lead. Under 1.8 cm. lead there 
some disagreement with the calculated pulse height distribution curve but this 
can attributed the production showers the lead. showers 
would obscure the presence small number particles unusually high 
ionizing power, such exist. 


Introduction 


recent years the measuring the ionizing power single cosmic ray 
particles has attracted careful attention. explaining some unpublished re- 
sults Weisz and Ramsay, Swann (11) has pointed out the importance 
statistical variations due the fact that appreciable portion the energy 
loss due ionization carried away secondaries which are the energy 
range thousands electron volts. chamber reasonable size the prob- 
ability occurrence these secondaries such that very poor statistics should 
expected the size ionization pulses due single particles. Alichanian, 
Alichanow, and Nikitin (1) and Nikitin (6) interpret results measurements 
the ionizing power particles indicating the existence positive and 
negative particles the soft component which are heavier than the meson, 
fact the latter suggests particles which can classified three mass groups 
considered worthwhile remeasure the ionizing power cosmic ray particles 
sea level whether these distinctions could confirmed whether 
ionization measurements would reveal anomalies which might indicate the 
existence appreciable number heavy mesons (3) the so-called 
particles reported recently different observers (7, 9). 


General Principles and Method Measurement 


The specific ionization strongly ionizing particles such a-particles 
slow protons usually measured ionization chamber filled with argon 
about atmosphere pressure. Owing the absence gas amplification 
such chambers, essential that the energy spent the particles within the 
vessel should sufficiently high. Otherwise, the resultant small ionization 
pulses would masked amplifier noise due Shot and Johnson effects and 
tube fluctuations. Ordinarily the collection about 5000 electrons the col- 

received April 18, 1951. 
Contribution from National Research Laboratories, Ottawa, Canada. 


Postdoctorate Fellowship under the National Research Council, Ottawa, 


| 
3 
i 
| 
j 
+ 


496 CANADIAN JOURNAL PHYSICS. VOL. 


lecting electrode ionization chamber low capacity and connected 
suitable linear amplifier gives pulse from which the energy dissipated the 
particle producing ion pairs may determined with reasonable accuracy. 
This involves expenditure energy 150 kev., since requires about 
ev. produce ion pair the gaseous volume. The relativistic particles 
cosmic radiation having specific ionization the order about ion 
pairs per cm. argon N.T.P. can therefore hardly detected non- 
multiplying chamber, unless the mean path length made inordinately long 
the gas pressure extremely high. The latter again leads other complica- 
tions revealed Swann’s experiment (10). 


The intrinsic merit proportional counter becomes evident its ability 
measure very low ionization. Filled with suitable mixture gases, mostly 
argon and methane, proportional counter can operated voltage which 
gives charge multiplication is, however, generally maintained 
the multiplication range between 10? and which characterized 
strict proportionality between the primary ionization and the collected charge. 
Thus, the formation electrons even less within the proportional counter 
responsible for collection approximately 5000 electrons the input grid 
linear amplifier, resulting output pulse whose amplitude measur- 
able with reasonable accuracy. 


the realm cosmic rays the proportional counter was first utilized 
Korff (4) when attempted measure the number highly ionizing par- 
ticles the upper layers the atmosphere. investigations, however, 
the cosmic ray particles crossed the proportional counter all directions, 
traversing different path lengths and producing pulse height distribution 
which did not characterize the ionization power individual particles. 
therefore necessary use the proportional counter telescope system, 
means which narrow beam cosmic rays selected Geiger counters 
passed through proportional counter, the passage the particle and the 
number ion pairs formed individual particle being recorded simul- 
taneously. This method proportional counter telescope arrangement was 
adopted both Weisz and Ramsay and Alichanian a/. for the study 
the relative specific ionization the cosmic ray particles. only the 
method recording the ionization pulses from the proportional counter that 
their experimental arrangements differ essential details. the former ex- 
periment, the pulse produced each discharge was recorded photographically 
giving differential curve pulse distribution. the latter experiment the 
selection the pulses was performed means so-called biased multi- 
vibrator single amplitude discriminator enabling one measure the 
number pulses with the height exceeding certain value and thus giving 
integral curve the pulse distribution. 


Having obtained the size frequency distribution the ionization pulses due 
cosmic ray particles, the next step consists the calibration the system. 
Since the region strict proportionality proportional counter function 


. 
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initial ionization, a-particles can hardly used for standardizing the ion- 
izing power cosmic ray particles. more appropriate use for comparison 
cosmic ray mesons filtered through few centimeters lead which have 
ionization about and assign them unit ionization power. 


The number ion pairs created fast particle the counter detet- 
mined its specific ionization power and the path length traversed the 
particle within the gas the counter. When fast particle loses energy 
ionizing gas multiple process takes place: 


(a) primary parent ion pairs, usually consisting electron and positive 
ions, are formed; 


secondary ion pairs are created the electrons positive ions the 
primary pairs when such primaries have sufficient energy. turn these secon- 
daries may create additional secondary pairs. Williams (12) has given 
expression for the number ion pairs per cm. viz., 


where known the probable ionization, and are constants for given 
kind particle, and v/c, being the velocity the particle and the 
velocity light. 


The value the most probable ionization, lies between the values the 
primary ionization and the total ionization, which includes the secondary 
ionization produced electrons ejected from the parent ions. Generally the 
total number ion pairs created fast particle will fluctuate appreciably 
about the value This mainly due the statistical fluctuation the 
number these secondary electrons, which have appreciable ionizing power 
themselves. 


Landau (5) has calculated the size distribution fluctuations ionization 
produced fast particles traversing through thin layer matter. His 
curves show maximum corresponding the most probable ionization but 
indicate that considerable fluctuation the ionization produced any par- 
ticular particle expected. 


the present experiment the calibration was carried out according the 
method Alichanian al. viz., obtaining the fluctuation curve and deter- 
mining from the size the pulse representing the most probable ionization. 


Each experiment consists essentially four steps: 

(1) Plotting the fluctuation the ‘‘bias curve” for mesons particular 
arrangement the counter telescope after inserting cm. lead plate between 
‘the proportional counter and the bottom tray Geiger counters. 

(2) Plotting the total ionization curve” without the lead filter. 

(3) Decomposition the ionization spectrum the soft component 
subtracting the first curve from the second. 
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(4) Adjustment Landau’s fluctuation curve with the ionization spectrum 
the point downward slope the latter. 


analysis (4) expected furnish approximate value the ioni- 
zation power new particle, any. 


Experimental Arrangement 


The counter telescope system, schematically depicted Fig. consisted 
one proportional counter mounted vertically between two trays Geiger 
counters. Each tray contained three self-quenching Geiger counters the 


mm. 
TUNGSTEN WIRE 


CALIBRATION 
WINDOWS 


Fic. arrangement the proportional counter telescope. 


metal type, cm. long and 2.5 cm. diameter connected parallel one 
coincidence channel. The particles could enter the system within angular 
aperture about 20°, traversing nearly equal lengths within the proportional 
counter. The maximum difference path length was approximately 10%. 
Suitable provisions were made for inserting cm. lead plate between the 
proportional counter and the bottom tray and also for placing lead plates 
appropriate thickness above the top tray. 
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The arrangement the proportional counter was somewhat similar that 
adopted Alichanian al. The central wire was supported side tubes. 
This design more convenient use telescope system than the tradi- 
tional type with the central wire sticking out both ends. The cathode was 
copper cylinder cm. long and 8.5 cm. diameter. tungsten wire, 0.1 mm. 
diameter, served the anode and was held tension spiral spring 
each end. The cylindrical envelope was made pyrex glass, the tungsten leads 
being taken out through nonex glass sleeves. The re-entrant glass barrel was 
provided with two capillary side tubes which were closed means thin 
glass windows having stopping power equivalent about one centimeter 
air. allowed narrow pencil a-particles from polonium source pass 
into the counter either the vicinity the anode wire near the cathode 
cylinder. The inside the counter was treated with concentrated acetic acid 
for few hours, rinsed several times with distilled water, and finally washed 
with absolute alcohol. The counter was next connected vacuum manifold 
and degassed for few hours temperature 350°C. under continuous 
pumping. Finally, the counter was filled with pure mixture containing 25% 
argon and 75% methane pressure cm. mercury. The outer surface 
the glass envelope, excepting the lead-carrying sleeves, was coated with 
thick conducting layer colloidal graphite and grounded. Before assembling 
the apparatus, the pulse-formation characteristics the proportional counter 
were tested firing a-particles through the calibration windows and exam- 
ining the shape and duration the pulses means cathode ray oscil- 
lograph. 


The main requirement the electronic equipment determine coinci- 
dences between Geiger counters and the proportional counter function 
the height the pulses from the latter. Since the pulses due the proportional 
counter are usually about few thousandths the amplitude those due 
Geiger counters, necessary have additional stages amplification pre- 
ceding the mixer stage. The anode wire the proportional counter was directly 
connected the grid the first tube the linear amplifier, which consisted 
two stages with amplification factor about 2000. The amplifier was 
the resistance capacity coupled type, having frequency characteristic such 
that accepted and preferentially amplified the pulses the desired type. 


The selection the pulses according their amplitude was performed 
means biased multivibrator, which accepts only pulses larger than par- 
ticular size. The negative pulse from this selector stage was fed Rossi tube 
which was connected twofold coincidence with second Rossi tube, the 
latter being operated the coincidence pulse from the Geiger counter trays 
after the usual phase inversion. 


The Geiger counters each tray were connected parallel, since the total 
capacity the wire system and counting rate were small. The anode wires 
the counters were directly joined the grids the corresponding Rossi tubes, 
while appropriate negative potential was applied the cathode. This made 
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the resolving time very small, ensuring short duration the pulses and avoid- 
ance accidentals. However, since the reaction time mechanical recorder 
requires that the recorded pulse prolonged time the order about 
one hundredth second, neon tube coupled multivibrator delay circuit 
was used between the Rossi and output stages for directly recording the coin- 
cident counter pulses. The chief advantages such circuit are that every 
pulse reaching the grid the output tube equal amplitude and duration 
and operates extremely small power without having any extinguishing 
troubles associated with thyratron. The double and triple coincidence pulses 
were recorded separately pen recorder. 


Separate power packs were used for the linear amplifier and coincidence 
units. The negative bias the multivibrator was obtained from calibrated 
potentiometer across 45-volt Burgess super service battery. The Geiger 
and proportional counters were supplied from independent high voltage 
sources; the high voltage power source for Geiger counters was obtained from 
neon stabilivolt, while that the proportional counter was obtained from 
separate battery 45-volt units series and contained completely 
shielded box. 


feeding artificial pulses from signal generator, was possible deter- 
mine whether the linear amplifier and pulse analyzer had drifted during cer- 
tain course measurements. This calibration unit consisted completely 
shielded battery circuit, operated with microswitch and producing pulses 
various sizes which could selected will. Together with the calibration 
corresponding proportional counter pulse, artificial pulse equal 
amplitude the Geiger counter pulses could taken off from suitable point 
the circuit and coupled through condenser the grids the respective 
Rossi tubes produce exactly the same effect actual threefold 
coincidence through the Geiger counters and the proportional counter. 
careful check was made means cathode ray oscillograph sure 
that the calibration pulses were identical shape and duration those the 
proportional counter and the Geiger counters. The calibration unit was also 
instrumental providing accurate information about the multiplication factor 
the proportional counter and the sensitivity the multivibrator during 
particular series measurements. 


Experimental Results 


The experiments were carried out the National Research Laboratories 
Ottawa, which has altitude approximately 170 ft. above sea level. The 
counter telescope was installed small cabin with thin roof wood and 
asphalt. 


The main investigation can subdivided into three sets experiments. 
the first set attempt was made detect highly ionizing particles the 
soft component the cosmic radiation. The results are represented Fig. 
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Curve the fluctuation curve the relativistic mesons, was first obtained. 
For this set measurements the soft component was absorbed lead 
filter placed between the proportional counter and the bottom tray Geiger 
counters, shown sketch the inset. When the dial reading the cali- 
brated potentiometer across the bias battery below 2.5 divisions, i.e. with 
the bias the multivibrator less than —11.25 volts, the curve goes parallel 
the abcissa. With further increase negative bias (greater pulse-height 
discrimination) the ordinates show rapid decrease followed more gradual 
slope. The most probable ionization lies the dial reading 2.7 divisions 
the potentiometer. Since only relativistic mesons were being recorded this 
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Fic. spectra the hard and soft components cosmic radiation sea level. 


arrangement, the dial reading 2.7 divisions may equalled unit ioni- 
zation power along the abcissa scale. Next the arrangement was adopted 
which gave the ionization spectra constituting the total intensity cosmic 
radiation. The ionization spectrum the soft component was obtained 
subtracting curve from curve and represented curve comparison 
with the fluctuation curve calculated according theory and fitted 
with the point its downward slope shows that the experimental and 
theoretical curves coincide with each other. 


the second set experiments 1.8 cm. lead plate was placed above the 
top tray Geiger counters. This about the optimum thickness for the for- 
mation cascade showers and would represent first maximum the 
number particles observed. The results are shown Fig. explained 
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above, and represent the fluctuation curves for relativistic mesons and 
the total intensity under these experimental conditions. The difference curve 
however, this case does not coincide with the calculated Landau fluc- 
tuation curve, the latter being represented the broken line. The difference 
between and the Landau curve small but extends large values the 
abscissa scale. contrast the agreement found Fig. and Fig. the dis- 
agreement the results with this experimental arrangement seems have 
high probability being real even though individual points the curves 
not show high statistical accuracy compared the difference between and 


the Landau curve. 


The third set experiments was carried out with cm. lead plate above 
the top tray Geiger counters. This arrangement was chosen being most 
favorable for the study the ionization power penetrating particles which 
are emitted under conditions characterized Rossi’s second maximum. Inci- 
dentally such experiment might also clarify the question whether many 
mesons are generated simultaneously penetrating shower such narrow 
cone that they cannot detected the usual shower-detecting arrangement 
Geiger counters. The results are represented Fig. the previous 
cases, curves and represent the fluctuation curve relativistic mesons 
and the ionization spectrum the total intensity under these experimental 
conditions. Curve the ionization spectrum the soft component obtained 
‘by subtracting curve from curve The differential curve also coincides 
with Landau’s fluctuation curve fitted with the point its downward 


slope. 
Discussion Results 


The results the first set experiments show that the soft component 
cosmic radiation sea level consists mostly relativistic electrons. There 
indication the presence highly ionizing particles the type discovered 
Alichanian al. altitude 3250 meters. Rossi and Greisen (8) at- 
tempted explain the total intensity the soft component under the assump- 
tion that all electrons originate from meson decay and knock-on processes. 
must emphasized, however, that, account the finite thickness the 
walls the counters, the particles the low energy end the spectrum and 
those with high specific ionizations (small ranges) are automatically excluded 
from detection our apparatus. would interesting repeat the experi- 


ment with thin-walled counters. 


The results the second set measurements are interesting. Accepting 
real the difference between curve Fig. and the Landau curve (broken 
curve Fig. 3), these results indicate the presence either appreciable 
number showers, particles high ionizing power. The former seems 
most likely since the experimental conditions are appropriate and Dobrotin (2) 
has attempted explain similar deviation his distribution curve high 
altitudes being due large number air showers. few particles with 
ionizing power corresponding charge greater than would cause disagree- 
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ment with the Landau curve large values the ionizing power, but showers 
several particles and particles near the end their range would have the 
same effect. These results, therefore, are inadequate show the existence 
particles unusual ionizing power the cosmic ray flux ground level. 
However, they demonstrate the way which the statistical variation 
ionizing power and the presence showers make identification ionization 
measurements very difficult. 


The result the third set experiments shows that even under cm. 
lead there appreciable number electrons, created the 
knock-on process mesons. The differential curve shows that these electrons 
are predominantly relativistic. this set measurements there hardly 
any evidence the production pairs mesons with such small angle 
separation that they would escape detection Geiger counter systems. 
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DISTRIBUTION ZENITHALE RAYONNEMENT COSMIQUE! 


Par 


Sommaire 


Les mesures ont porté sur composante pénétrante qui est surtout composée 
mésons niveau sol. Deux bandes contigues ont été sélectionnées dans 
spectre différentiel; les impulsions des particules enregistrées 
300 410 Mev./c 410 510 Mev./c pour des mésons. distribution 
leur intensité été étudiée pour les angles zénith compris entre 180°. Les 
résultats indiquent que loi distribution zénithale des particules 
étudiées est différente celle qui exprime distribution zénithale spectre 


intégral. 


Introduction 


zénith été étudiée dés 1932 (1). est bien connu que, niveau mer, 
maximum rayonnement provient direction verticale que, 
par suite atmosphérique radioactivité des mésons, cette 
intensité decroit pour des directions éloignées zénith. loi empirique qui 
exprime distribution spectre intégral, toutes les particules 
qui subsistent aprés avoir été filtrées par écran absorbant, peut-étre repré- 
rayonnement observé dans une direction faisant angle avec zénith 
rayonnement prend une valeur peu supérieure (environ 2.2) dans cas 
des particules pénétrantes (4, 12) filtrées par une dizaine cm. plomb. 
Par contre, nous avons montrer, dans une expérience précédente (10) que 
cette loi s’applique pas dans cas spectre différentiel des mésons. 


des expériences qui vont étre décrites est préciser loi qui gou- 
verne distribution zénithale des mésons dont est comprise entre des 
limites déterminées. effet, est important connaitre distribution 
zénithale ainsi que spectre différentiel des mésons dans les problémes qui 
rapportent spectre primaire rayonnement cosmique, production des 
spectre différentiel des mésons qui est encore assez mal connu été étudié 

_expérimentalement collaboration avec Monsieur Rose, les résultats 
dernier travail seront publiés séparément. nous avons essayé d’é- 
tendre les mesures zénithales aux angles suivants: 90°, 120° 180° avec 
zénith, afin d’obtenir distribution des particules pénétrantes. 
effet, Powell ses collaborateurs (3) ont observé une composante rétrograde 
‘dans des émulsions photographiques exposées Jungfrau; par ailleurs Ritson 
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sol. serait intéressant d’étudier distribution angulaire ces 
afin d’aider sans doute déterminer L’expérience que nous avons 
effectuée permet pas confirmer cette composante rétro- 
grade dont est faible niveau sol. 


Dispositif Experimental 

Principe Méthode 

Nous avons une méthode d’absorption différentielle qui consiste 
enregistrer les particules qui peuvent traverser une certaine épaisseur plomb 
mais qui sont arrétées dans écran absorbant additionel. cette maniére 
sélectionne les particules dont parcours est ainsi déterminé: 
étudiée est déterminée par différentiel son emplace- 
ment dans spectre par Les particules qui sont ac- 
compagnées d’une gerbe sont éliminées fagon n’enregistrer que les par- 
ticules pénétrantes qui sont, pour plus grande part, des mésons parmi 
lesquels trouve probablement mélangé faible pourcentage protons 
dans les conditions fréquence des particules ainsi sélect- 
ionées est mesurée pour différents angles zénith loi expérimentale qui 
exprime distribution peut-étre établie. 


Description 

télescope double compteurs A.B.C triple coincidences (Fig. 
définit faisceau particules qui ont traversé les écrans 
tivement les particules qui ont traversé les écrans chacun épais- 
seur 7.5 cm. plomb. 


Les groupes compteurs enregistrent les gerbes. Les compteurs les 
écrans plomb sont montés sur solide fer qui permet 
télescope selon angle quelconque avec zénith. 


triple coincidence A.B.C constitue principale d’un circuit 
groupe maitre. L’impulsion maitresse, provenant télescope, est elle-méme 
coincidence avec chacun des groupes compteurs Les coin- 
cidences suivantes sont alors enregistrées: A.B.C, A.B.C.D, A.B.C.E 
A.B.C.G desquelles peut-étre déduit les anti-coincidences 
qui correspondent aux particules qui sont pas accompagnées d’une gerbe 
qui s’arrétent dans ainsi que les anti-coincidences A.B.C.D-(E.G) 
qui correspondent également aux particules solitaires mais qui, cette fois, 
dans Deux bandes d’énergie contigues sont ainsi dé- 
coupées dans spectre différentiel enregistrées simultanément fonction 
zenith. 


geométrie des compteurs été étudiée répondre, autant que 
possible, aux conditions suivantes: 
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(a) L’ouverture solide soutenu par télescope doit aussi 
faible que possible, bien que compatible avec une fréquence comptage 
raisonnable. Une finesse relative solide défini par télescope réduit 
mesurer distribution zénithale des directions voisines 


Fic. Disposition des compteurs des écrans absorbants. 


zontale. d’ouverture minimum télescope mesure 20° 
latéral 40° qui autorisé des mesures pour les angles zénith 80° 
sans que télescope soit traversé fois par les particules venant deux 
directions azimuthales opposées. 


(b) Les compteurs anti-gerbes ont été disposés vue bonne efficacité 
tant pour les gerbes produites dans les écrans supérieures que pour les gerbes 
Les deux groupes compteurs sont adjacents faisceau 
rayonnement déterminé par télescope, qui les rend sensibles méme aux 
gerbes dont les particles sont peu divergentes. L’ecran d’une épaisseur 
2.5 cm. plomb accroit probabilité production gerbes par les parti- 
cules gerbigénes arréte les gerbes les plus molles. Nous espérons 
obtenir ainsi une discrimination bonne entre les composantes 
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molles pénétrantes rayonnement cosmique (9) dans domaine des éner- 
gies étudiées. 

Afin d’avoir une protection suffisante contre les gerbes latérales, pour les- 
quelles fréquence relative est importante aux grands angles, les compteurs 
sont placées maniére couvrir complétement les compteurs Les écrans 
augmentent probabilité d’enregistrer les électrons les photons 
accroissant densité des gerbes dessus des compteurs Des circuits sup- 
plémentaires d’amplification réduisent une valeur négligeable 
des compteurs qui sont paralléle. méme dispositif d’amplification est 
utilisé pour chacun des groupes dans lesquels les compteurs 
sont aussi connectés 


(c) plusieurs avantages enregistrer simultanément deux bandes 
spectrales contingues. particulier, une estimation expérimentale 
diffusion dans différentiel peut effectuée. Mais 
deux bandes contigues est surtout intéressante dans détermination 
spectre différentiel pour lequel nous obtenons ainsi recouvrement des deux 
bandes contigues mesurées. 


part, peut étre exercé permanent sur des 
compteurs car nombre des particules manquées par ces compteurs, mais 
qui sont enregistrées par les compteurs doit rester trés faible dans les con- 


ditions normales. 


Discussion 


Les différentes causes qui peuvent d’erreurs les résultats dont 
doit étre tenu compte seront discutées dans les lignes qui suivent: 


(a) Fond continu 

l’on retire les écrans différentiels une certaine fré- 
gine fond continu plusieurs causes parmi lesquelles les plus impor- 
tantes sont probablement les suivantes: 


(1) Inefficacité des groupes compteurs une particule traverse 
compteur dans une région faible efficacité, pendant 
temps mort, c’est dire moins quelques sec. aprés passage d’une 
particule précédente, donnée par compteur est faible ampli- 
tude. remédier cet inconvénient, une amplification convenable 
des impulsions est appliquée méme maniére que pour les compteurs anti- 
gerbes. Toutefois, cette cause d’inefficacité des compteurs pas seule, 
peut aussi produire des impulsions retard etc. subsiste toujours 
fond continu aux compteurs. 


(2) Absorption dans les parois des compteurs.—Cette absorption est calculable 
représente environ 10% fond continu total. D’ailleurs, elle n’est pas 
effet parasite provoque véritables anti-coincidences. 
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(3) Coincidences fortuites triples A.B.C.—Le pouvoir résolution 
est: sec. nous écrivons que fréquence des triples 
coincidences fortuites est premiére approximation: 
étant fréquence d’un double compteur télescope, nous trouvons qu’il 
produit coincidences fortuites par heure. Cette valeur est négligeable 
comparaison avec fréquence fond continu qui est 0.5 
anti-coincidences par heure. 


(4) Inefficacité des compteurs gerbe latérale faible exten- 
sion peut déclencher les compteurs sans atteindre les groupes 
Une gerbe type elle était manquée par les compteurs serait enregistrée 
comme anti-coincidence par conséquent serait indiscernable des particules 
pénétrantes arrétées dans les écrans différentiels. Toutefois, nous avons 
que dispositif anti-gerbes est certainement convenable. Une 
peut étre faite. Dans les conditions les moins favorables, quand 
télescope est incliné angle 80° avec zénith, fréquence des gerbes 
latérales est maximum taux comptage télescope est son minimum. 
Dans ces conditions, télescope est déclenché fréquence 7.1 coincidences 
par heure dont 2.0 coincidences sont accompagnées gerbes. taux enregistré 
pour fond continu est 0.13 coincidence par heure. nous supposons que tout 
fond continu est aux gerbes latérales nous trouvons que sur deux gerbes 
par heure est manqué seulement 0.13, c’est dire que 
groupe anti-gerbes est bien supérieure 94%. 
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Généralement, dans les mesures analogues celles que nous avons effec- 
tuées, fond continu est mesuré pour direction verticale pourcentage 
par rapport totale est supposé d’une valeur constante pour tous 
les angles. pourcentage est utilisé pour corriger fond continu les anti- 
coincidences enregistrées quand les écrans différentiels sont place. Toute- 
fois, une étude systématique fond continu pour differents angles zénith 
indique que rapport étant fréquence des anti-coincidences 

m 

celle des mésons traversant télescope, s’accroit (voir Fig. 
quand verticale. Cet accroissement relatif fond continu 
pour les grands angles, particulier pour 75° 80°, pourrait produit 
non seulement par plus faible intensité J,, mais aussi par présence prés 
composante trés molle absorbée dans les parois des compteurs. 
Cette derniére supposition est d’ailleurs appuyée par fait que nombre des 
anti-coincidences enregistré par groupe est plus faible que celui indiqué 
par groupe pour les angles 75° 80° avec zénith. 


Dans les résultats des mesures qui sont donnés plus bas, fond continu, 
mesuré pour chaque angle, été soustrait fréquence par heure des 
particules pénétrantes arrétées par les écrans différentiels 


(b) Diffusion 

Dans supérieur diffusion est négligeable car nombre 
des particules qui sont écartées faisceau défini par télescope est, moins 
partiellement, compensé par nombre des particules qui originellement 
étaient dehors solide défini par télescope qui, par suite 
diffusion, sont orientées solide. Dans les écrans 
méme effet peut produire, plus, ces écrans étant faible épaisseur, 
diffusion des protons des mésons rapides n’y est pas considérable. Enfin, 
les mesures qui sont décrites ici n’ayant pas pour objet détermination d’une 
intensité absolue, diffusion été négligé, sans inconvénient croyons- 
nous, pour les absorbants P3. 


Dans les écrans des particules peuvent aussi diffusées dehors 
faisceau mais, cette fois, est préjudiciable car, telles particules 
sont fortement déviées, elles peuvent passer dehors des compteurs 
aussi confondues avec les particules absorbées dans écrans différ- 
entiels. Afin parer autant que possible ceci, les groupes compteurs 
couvrent largement solide défini par télescope. 


L’étude géométrie des compteurs montre que, dans 
cas typiquement défavorable, les particules diffusées devraient avoir 
angle 30° moins avec leur direction primitive sortie pour 
étre manquées par les compteurs est improbable qu’un événement 
genre produise avec une fréquence suffisante pour altérer résultat des 
mesures. Mais n’en pas méme qui concerne diffusion dans 
second écran différentiel. Dans cas, méme une faible déviation des parti- 
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cules dans peut les écarter surface utile des compteurs 
plus, une particule déja déviée par premier écran peut encore étre diffusée 
dans second. L’estimation proportion des particules diffusées qui 
manque groupe peut étre fait comparant les indications des groupes 
Afin d’éviter une hypothése sur courbe spectre différentiel, 
méme bande spectrale est pour les deux groupes compteurs, ceci est 
obtenu changeant Les valeurs enregistrées pour les fré- 
quences sont respectivement corrigées fond continu, correction 
qui inclus celles attribuables aux inefficacités des groupes compteurs 
différence qui subsiste entre AJp qui mesurent méme bande 
spectrale, peut étre attribuée diffusion. Nous trouvons ainsi que 
valeur est environ plus faible que celle trouvée pour Dans nos 
résultats valeur été corrigée pour cet effet. 


(c) Gerbes 

Les gerbes sont enregistrées par les compteurs qui sont coincidence avec 
télescope. pouvoir pénétration des gerbes peut étre estimé par 
nombre celles qui traversent les écrans Quand télescope 
est horizontal, fréquence des gerbes enregistrées indique limite supérieure 
fréquence des gerbes latérales, car c’est dans ces conditions que 
enregistre maximum des gerbes qui sont principalement direction 
verticale. Les courbes représentées sur figure montrent fréquence des 
gerbes fonction zénith. peut remarquer que rapport des 
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gerbes pénétrantes sur nombre total des gerbes est élevé. Par exemple, pour 
direction verticale trouve qu’il environ 55% gerbes pénétrantes, 
l’on déduit les gerbes latérales. plupart ces gerbes pénétrantes sont 
probablement provoquées par les électrons collision des mésons. Elles repré- 
sentent 4.3% composante verticale des mésons, valeur qui 
est méme ordre que celle trouvée d’autres méthodes (2-13). Pour 
les grands angles avec zénith fréquence relative des gerbes pénétrantes 
décroit nombre des gerbes latérales devient prépondérant. 


(d) Correction barométrique 

valeur correction barométrique effectuer sur les mesures est assez 
bien connue qui concerne spectre intégral, n’en n’est pas méme pour 
spectre différentie]. travail récent (11) méme montré que parfois cor- 
rection pourrait étre positive, tant6t négative suivant domaine 
d’énergie considéré. Les mesures qui sont décrites ici ont été effectuées durant 
une période plusieurs mois pendant laquelle pression atmosphérique est 
demeurée relativement stable. moyenne pression atmosphérique 
pour chacune des mesures s’établit comme suit: 

30° 60° 75° 80° 

millibars 1007.7 1008.5 1007 1006 1009 

pression atmosphérique moyenne pendant durée des expériences été 
1007.6 millibars. Par suite faible déviation pression durant 
chacune des mesures autour pression moyenne (au maximum 1.4; 
1.6 millibars) par suite sur valeur coefficient baromé- 
trique qui concerne les mesures différentielles, n’a pas été effectué 
correction barométrique sur les résultats. 


Distribution zénithales pour les angles 90° 

Les mesures ont été effectuées Ottawa dans batiment construction 
légére pour les angles zénith suivant: 0°, 30°, 60°, 75°, 80°, 90°, 120° 180° 
dans direction Les résultats obtenus pour les angles 90°, 120° 
180° seront exposés dans paragraphe suivant. 


Durant ces mesures des écrans été maintenu constante égale 
nous supposons que toutes les particules sélectionnées sont des mésons, leurs 
impulsions s’etendent, d’aprés les courbes parcours fonction 
pulsion (7) 300 410 Mev./c 410 510 Mev./c. 


Les résultats sont résumés sur les tableaux II, qui indiquent nombre 
particules absorbées dans chacun des écrans différentiels 
respectivement, pour les angles zénith. Dans tableau III 


log cos 


varierait distribution était représentée par prenant comme 
référence verticale seulement angle avec énith. figure montre 
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comment les points expérimentaux placent par rapport courbe théorique 
Nous pouvons voir que prés verticale radiation décroit trés 
rapidement fonction zénith que cette décroissance est plus 
lente quand s’écarte direction verticale. Cet effet est particuliérement 
visible pour les particles absorbées dans premier écran différentiel, c’est 
dire pour bande d’énergie moins élevée. Ainsi que nous venons voir, 


TABLEAU 
PARTICULES QUI SONT ARRETEES PAR DIFFERENTIEL 


30° 60° 75° 80° 


TABLEAU 
PARTICULES QUI SONT ARRETEES PAR DIFFERENTIEL 


652 166 0.086 0.025 


TABLEAU III 


log 


30° 60° 75° 80° 
400 Mev./c) 


510 Mev./c) 


est impossible tracer une droite par les points courbe log 
fonction log cos parce que pas constant. Par contre, les points qui 
représentent log fonction log sin placent sur une ligne 
droite. Ainsi, nous obtenons loi empirique: 


sont des constantes ayant les valeurs suivantes: 


if 
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conclusion, les résultats montrent une forte décroissance prés zénith 
une tendance plus d’isotropie pour les directions éloignées verticale, 
particulier, dans cas des mésons 300 410 Mev./c. Les phénoménes 
qui pourraient produire cet effet d’isotropie sont diffusion des particules 
long leur trajectoire dans désintégration vol mésons 
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théorique les points expérimentaux correspondant aux mesures 


lourds donnant lieu des mésons plus légers les désintégrations nucléaires 
dans lesquelles des mésons sont créés. D’autre part, ces deux derniers phéno- 
pourraient produire une composante rétrograde provenant 
direction 90° avec zénith. 


Distribution pour les angles 90° 

Les mesures par méthode des anti-coincidences sont délicates dans cas 
des angles 90° car, dans ces derniéres mesures, télescope est aussi dé- 
clenché par les particules ayant des directions opposées celles des particules 
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étudier. Néanmoins, des mesures préliminaires ayant donné des résultats 
encourageants Pic Midi, France, essai été tenté niveau 
sol. Comme pénétration radiation rétrograde pas connue, 
seur absorbant été variée entre cm. plomb. Mais, 
afin conserver une certaine discrimination entre les particules penétrantes 
les particules gerbigénes, les écrans intermédiaires chacun 2.5 
cm. d’épaisseur, ont été maintenus leur place. Par des mesures croisées, dans 
lesquelles les deux écrans différentiels étaient retirés, nous avons essayé 
d’estimer flux ascendant. 


(1) Mesures 90°.—Le tableau donne les résultats des mesures effectuées 
sous cet angle zénith. représente fréquence des particules qui 
télescope sans étre accompagnées par une gerbe. Nous pouvons voir, d’aprés 
tableau, que des particules direction horizontale est peu affectée 


TABLEAU 


EFFET DES ECRANS ABSORBANTS SUR COMPOSANTE HORIZONTALE 90°) 


~ 
& 


COCO 
o 


nombre par heure. 


Rogozinski dans une étude composante horizontale. qui concerne 
n’est pas possible dire cette augmentation est due 
neutre, venant une particle chargée. 


(2) Mesures 120° 180°.—Ces deux séries mesure présentent pas 
réellement significatif avec sans les divers absorbants. Aussi, nous 
avons effectué une autre série mesures méthode plus sensible 
basée sur détection des mésons qui desintégrent dans absorbant 
donnant naissance électron. Cette méthode d’ailleurs été utilisée par 
Ritson (8) dans des mesures analogues. figure montre schéma 
télescope A.B est coincidence retardée 1.2 47.2 usec., 
avec les compteurs C.D eux-mémes double coincidence. Les coincidences 
fortuites fond continu sont abaissés une valeur raisonnable 
par des circuits d’anti-coincidence. Ainsi sont enregistrés les mésons qui aprés 
avoir traversé télescope désintégrent dans brique carbone. Cette 
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méthode présente d’éliminer toutes les particules qui sont pas 
des mésons bien des mésons donnant naissance des mésons provenant 
direction zénithale choisie. fond continu est déterminé 
retirant brique carbone. est, dans nos mesures, égal 0.42 0.10 


ECHELLE 


Fic. Disposition des compteurs carbone dans coincidence 
différée. 


coincidences par heure. trouve pour 180° une fréquence 0.45 0.06 
coincidences par heure quand brique carbone est place. Ainsi cette 
mesure n’est pas plus concluante que précédente. résumé, ces expériences 
montrent que flux ascendant, s’il existe, est faible 
niveau sol. 

Conclusions 


Les mesures que nous avons effectuées Ottawa sur distribution zénithale 
des particules pénétrantes parcours compris entre cm. 22.5 cm. 
entre 22.5 cm. plomb indiquent que ces particules suivent pas loi 
habituelle mais que distribution leur intensité peut étre repré- 
sentée par loi empirique: 


dans laquelle est environ égal prend les valeurs suivantes: 
1.47 pour les particules les moins énergiques 1.61 pour seconde 
bande étudiée. fait plus intéressant, mis évidence par les mesures, est 
d’une décroissance rapide prés zénith une 
tendance plus d’isotropie prés bande spectrale plus 
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faible énergie présente ces caractéristiques d’une maniére beaucoup plus mar- 
quée que seconde bande, notamment pour les angles zénith 75° 80°. 
peut d’ailleurs tel résultat s’il une production locale 
particules pénétrantes, bien des désintégrations vol mésons masse 
supérieure méson bien encore des particles trés déviées leurs trajec- 
toires primitives par diffusion travers 


qui concerne radiation rétrograde, son intensité est trés faible 
niveau mer des expériences complémentaires altitude seraient 
désirables afin d’en déterminer 
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THE SOLUTION X-RAY ACTIVATION CURVES FOR 
PHOTONUCLEAR CROSS SECTIONS! 


Abstract 


method presented for the computation photonuclear cross sections 
from their X-ray activation yield curves. based the photon differences 
between successive Schiff representations bremsstrahlung spectra and assumes 
that the activation curves have smooth first and second derivatives. Artificially 
prepared activation curves have been satisfactorily solved this method, and 
the published photonuclear activation curves determined this laboratory have 
been reanalyzed. The resulting cross section curves are essential agreement 
with those originally determined; however the new values are believed 
more reliable. New constants for two proposed relationships between 
peak energies and atomic mass number have been determined. 
Appended the paper are tables the Schiff bremsstrahlung spectra for 
maximum photon energies between and Mev. and also special functions 
for cross section computations between the same energies. 


Introduction 


1948, Baldwin and Klaiber (2) published analysis their activation 
curves for the reactions and The cross section curves 
computed exhibited peaks with half-widths few Mev. Goldhaber and 
Teller (7) and Levinger and Bethe (19) then proposed theories account for 
certain characteristics such peaks the result resonance dipole absorp- 
tion photons nucleus, followed de-excitation through various com- 
peting modes decay. The predominant mode should neutron emission. 
Interest the shapes the cross section curves for photonuclear reactions 
was greatly stimulated these theories. 


For such investigations the ideal instrument would generator mono- 
chromatic y-rays variable energy. Many and reactions have been 
studied with y-rays from the reaction (12, 33). Little can 
concluded about the cross section shapes from such studies, since these y-rays 
are mixture 14.5 and 17.5 Mev. photons and are produced from narrow 
resonances the proton bombarding energy. Monochromatic y-rays energy 
greater than Mev. are produced from the reaction H*(p, The proton 
energy resonance for this reaction much broader, that one could hope 
use for construction such generator. However, this has not yet been 
done for the study photonuclear reactions. 


The only other source high energy y-rays bremsstrahlung produced 
when electrons strike target after being accelerated betatron synchro- 
tron. This process gives continuous spectrum photon energies. difficult 
sort out the photonuclear effects produced photons different energies 
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from given irradiation such X-ray beam. principle, the energies 
all the disintegration products for each large number photonuclear 
events could observed, then frequency distribution these events 
function photon energy could plotted. approximate knowledge the 
shape the bremsstrahlung spectrum would then suffice for the determination 
the cross section curve. far has been possible apply this procedure 
only such reactions 4a), where all the reaction products are a-par- 
ticles which can observed nuclear emulsion (25). 


Usually, however, possible measure only the yield one the dis- 
integration products reaction (such neutrons, protons, residual radio- 
active nuclei) function the dose and maximum energy the brems- 
strahlung spectrum which produced them. The cross section curve for the 
reaction can determined this case only both the photonuclear yield 
and the shape the bremsstrahlung spectrum are known functions the 
maximum energy. 


Johns, Katz, Douglas, and Haslam (13) have published curves based the 
Schiff bremsstrahlung spectra (1, 30) and outlined method for the compu- 
tation cross section curves from the corresponding X-ray activation curves. 
This method will referred hereafter the ‘‘total method. 
open the objection that progressive smoothing the cross section curve 
must carried out the course computation, and this some extent 
subjective and hence arbitrary process. The present authors have de- 
veloped method treatment, suggested Penfold, which com- 
putationally simpler and less arbitrary than the total spectrum method. This 
will designated the ‘‘photon method. 


II. Computational Functions 


The intensity the bremsstrahlung spectrum the forward direction 
according Schiff (1, 30) proportional 
where 
and 


these expressions the rest energy the electron, the photon energy, 
k/W, the total energy the electron, and the atomic number 
the target material. 


Schiff originally proposed (1, 30) that the constant should have the value 
191. However, from recent reanalysis bremsstrahlung theory (31), has 
obtained the value 111. Equation (1) represents the shape the brems- 
strahlung spectrum the forward direction only for infinitely thin targets, 
and the relationship thus cannot used stands give the bremsstrahlung 
spectra obtained from betatrons and synchrotrons, which have thick targets. 
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When Schiff’s equation, its angular dependent form (31), integrated 
over all angles, the resulting total bremsstrahlung curve practically indis- 
tinguishable shape from that Rossi and Greisen (29) and Heitler (11). 
Equation (1), with 191, actually compromise shape between the 
total bremsstrahlung spectrum and that the forward direction calculated 
with 111. Such compromise represents better the actual bremsstrahlung 
spectrum the forward direction from thick targets, where multiple traversals 
and scattering electrons give additional contribution the lower energy 
portion the bremsstrahlung curve. 


Koch and Carter (16) have measured the distribution bremsstrahlung 
photon energies determining the energies electron pairs produced from 
19.5 Mev. rays cloud chamber. Their resulting intensity spectrum 
reproduced Fig. The old and new Schiff curves computed from Equation 
(1) have been fitted the experimental points have the same number 
photons between and 19.5 Mev., where the results Koch and Carter 
are most reliable. may seen that the experimental data lie considerably 
above the curve 111 lower photon energies and that the curve 191 
gives somewhat better fit the data. similar experiment using 322 Mev. 
rays, Powell, Hartsough, and Hill (28) obtained excellent agreement with 
bremsstrahlung shape calculated Christian and modified account for the 
effects their target thickness. The modified Christian bremsstrahlung shape 
nearly identical with that calculated from Equation (1) using 191. 


Johns al. (13) have investigated the changes introduced into the deter- 
mination the cross section curve the total spectrum 


INTENSITY 


200 


RELATIVE 


PHOTON ENERGY (MEV) 
Fic. experimental results Koch and Carter the shape the bremsstrahlung inten- 


sity curve. old and new theoretical curves 191 and 111) have been normalized 
include the same number photons between and 19.5 Mev. the experiment. 
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method when radically different types bremsstrahlung curve were used, 
including extreme cases far outside the experimental uncertainties involved 
the use Equation (1) with 191. The relative shapes the cross 
section curves are quite different, but the peak positions are nearly 
the same and the integrated cross sections are practically unchanged. 


have therefore assumed that Equation (1), with 191, gives ade- 
quate representation the shape the bremsstrahlung spectrum obtained 
from betatrons and synchrotrons. Experimental results check this shape suffi- 
ciently well that photonuclear cross section curves, computed the photon 
difference method, should changed only slightly when accurate experimental 
shapes become available. the following calculations the value (for 
platinum target) was used. evident from Equation (1) that large changes 
atomic number correspond small changes that the tables prepared 
for this paper should applicable any heavy target. 


The “standard” irradiation conditions adopted this laboratory require 
that sample placed the X-ray beam behind cm. Lucite and that 
X-ray dosages measured Victoreen ionization chamber placed the 
same position. This thickness plastic material necessary order bring 
the bremsstrahlung photons into equilibrium with their secondary electrons. 
should noted that this amounts arbitrary definition the ‘‘roent- 
asa unit X-ray dose. different from the “roentgen” defined through 
the use lead converting material (24), since many more electron pairs 
are produced the lead than the light elements Lucite. Mayneord (21) 
has defined the the extrapolation zero thickness the linear 
portion the depth dose curve Lucite. Such about 19% 
smaller than the cm. roentgen for Mev. rays. 


This method monitoring X-ray dosages has been compared the Lax 
chamber method and others McElhinney and Siewers (23), using brems- 
strahlung energies Mev. Very good relative agreement was obtained 
for X-ray energies Mev., the region which the secondary electrons 
are brought into full equilibrium with the primary beam. 


thus necessary modify the Schiff curves take account the photon 
absorption the donut walls and monitor the betatron and for the cm. 
Lucite. They must then multiplied the response curve the Victoreen 
r-meter. This was done Johns (13), and the resulting spectrum curves 
were normalized that each corresponded 100 roentgens. These curves 
will designated P(E, Eo), where the photon energy (Mev.), and 
the betatron operating energy corresponding the upper limit the 
spectrum. 


The saturated activity that would induced sample after infinitely 
long irradiation the X-ray beam 


Eo 
0 
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where the saturated activity per mole parent the cross section 
for the reaction barns; and P(E, Eo) the number photons per cm.? per 
interval per unit time. The numerical factor Avogadro’s number divided 
convert the cross section units from cm.? barns. The unit time 
may cancelled out the above expression, which case becomes the 
number reactions taking place per 100 per mole parent isotope. 
here defined gm-atom the isotope. 


Replacing the integral Equation (2) summation over series 
energy intervals and average cross sections width Mev., one obtains 
n=0 
a(Eo) = 0.6023 (n + P [Eo (n + de, Eo], (3) 
n = Eo/e 
where P(E, Eo) must now expressed units photons per per Mev. 
interval per 100 r., and must multiple The differences between 


such activities intervals Mev. are 


AP(E, Eo) P(E, Eo) P(E, 
solve for the cross sections this rewritten 
where AA(Eo) (6) 


Equation (5) suited for the rapid computation photonuclear cross sec- 
tions from their activation curves. Activation differences are taken and multi- 
plied the factors Equation (8) obtain the quantities The 
cross sections may then computed tabular form illustrated Table 


TABLE 


ILLUSTRATION OF THE COMPUTATIONAL PROCEDURE FOR THE CALCULATION OF A 
PHOTONUCLEAR CROSS SECTION FROM AN ACTIVATION CURVE 
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Under the series energies the quantities are entered. 
When the cross section energy determined, multiplied the 
values Eo) which are subtracted the appropriate columns the 
table. then the sum all such products under For accurate com- 
putations, account should taken the fact that the first activation differ- 
ence above the threshold the reaction does not correspond width Mev., 
and appropriate adjustments should made. 


order implement this procedure was necessary obtain tables 
the functions and Eo). was decided choose 0.5 Mev. The 
basic curves P(E, Eo) obtained Johns were cross-plotted give lines 
constant photon intensity versus which were smoothed. interpolation 
this cross-plot new set curves Eo), expressed units photons per cm.? 
per 0.5 Mev. interval per 100 r., were obtained for 0.5 Mev. The 
quantities AP(E, were then obtained subtraction. Because the brems- 
strahlung spectral shape changes slowly function Eo, was then insisted 
that AP(E, should smooth functions photon and betatron energies. 
This smoothness was achieved series cross-plots against the relevant 
variables. The new AP(E, Eo) were then corrected small smoothed factors 
that they could summed into good agreement with the original curves 
P(E, 

The resulting summation given Appendix this paper. The quan- 
tities defined Equation (8), and Eo), defined Equation (7), 
were computed from the smoothed AP(E, and are given Appendices 
and III, respectively. Only the second decimal place significant these 
tables, but the third has been retained for reasons smoothness. 


III. Derivatives the Activation Curves 


Individual activation points are usually determined only accuracy 
5%, owing errors the determination the initial activity the 
sample, the X-ray dose delivered the sample, and many cases the 
effective betatron operating energy.* therefore necessary smooth the 
activation curve. This can produce considerable distortion the computed 
cross-section curve the smoothing not done satisfactory manner. 
remove the arbitrary features this step necessary obtain criteria 
govern the smoothing procedure. 


The various methods analysis which have been applied experimental 
activation curves have all yielded peaked cross section curves (see for example 
13, 14, 15, 22, 26). peaks thus seem characteristic 
feature photonuclear cross section curves. would instructive examine 
the characteristics the derivatives the activation curve corresponding 
such peaked cross section. 


*The energy control circuit used the University Saskatchewan betatron effective Mev. 
this region the energy believed constant given setting +0.1 Mev. and calibrated 
+0.2 Mev. Above Mev. self-expansion the electron beam used which best holds constant 
Mev. with equal error calibration. 
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The following analysis indicates the general nature these derivatives. 
For mathematical simplicity cross section curve assumed the form 


photon spectrum, which when normalized constant 
dose (assuming the r-meter response linear), gives 


P(E, Eo) (E, (10) 
a 


Equation (11), the activation curve, plotted Fig. Note that although 
the cross section still quite large for energies above Mev., the activation 
curve has leveled off this region and insensitive the tail the cross 
section curve. 


The first and second derivatives the activation curve are then 
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Fic. activation curve and its derivatives calculated for assumed peaked cross section 
curve and photon spectrum. (a). The assumed cross section curve. 
The activation curve. (c) The first derivative. (d). The second derivative. 
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These are plotted Figs. 2(c) and 2(d). They are smooth functions, and the 
derivatives experimental activation curve resemble them gen- 


eral shape. 


felt that the best criterion fit experimental activation points 
curve having smooth first and second derivatives the general shape those 
Figs. 2(c) and 2(d). With the exception the reactions 
which will discussed later, has been possible this procedure fit 
activation curve the experimental points, determined this laboratory for 
several photonuclear reactions, within the experimental errors. practice 
good smooth curve was drawn through the activation points. The activation 
differences (which are proportional the first derivative) were read 
off this curve. These were then smoothed such way that their derivative 
was itself smooth and also such that the area under this smoothed curve was 
the same the area defined the activation difference points. 


The above procedure has the following pragmatic advantage. When experi- 
mental activation points were smoothed this manner, the progressive com- 
putation the cross section curve according Equation (5) yielded smooth 
succession cross section points. Any small oscillations these points were 
found damp out rapidly. The introduction during the smoothing procedure 
discontinuities extra inflections into the activation difference curve 
general caused much larger oscillations which took longer damp out. 
will shown later, small bulge this curve did not cause oscillations 
the calculation the cross section curve for the reactions 


The smoothing procedure has also been tested through the solution arti- 
ficial activation curves. number geometrical cross section curves were 
assumed one (L.K.) and the corresponding activation curves were 
computed the inverse the total spectrum method, using the curves 
published Johns (13). allowance made only for errors reading 
these curves, then the resulting activation points are subject uncertainties 
about two per cent. The other (A.G.W.C.) was unable predict 
inspection any distinctive features these cross section curves from their 
activation curves, apart from the obvious fact that they were peaked. 


Activation differences were taken from the artificial curves, and these were 
smoothed according the above procedure. The original cross section curves, 
together with the points then computed for these curves, are shown Figs. 
and may noted that the sharp bends the original curves were 
quite well reproduced, although the very sharp bend the curve Fig. was 
somewhat smoothed over. each case the initial rise the cross section curve 
was very well reproduced. Owing the insensitivity the activation curves 
the the cross section curves, these tails were not well duplicated. 


The authors are thus convinced the validity the smoothing procedure 
outlined above. 
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(ARBITRARY UNITS) 


SECTION 


CROSS 


10 “92 14 16 is 20 22 24 
PHOTON 


solid line assumed triangular cross section curve, which was used the 
computation activation curve from the spectra Johns al. The points were recomputed 


for the cross section curve the photon difference method. 


(ARBITRARY UNITS) 


SECTION 


ENERGY 
Fic. The solid line assumed quadrilateral cross section curve, which was used the 
computation activation curve from the spectra Johns al. The points were recomputed 
for the cross section curve the photon difference method. 
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(ARBITRARY UNITS) 


SECTION 


PHOTON 


Fic. solid line assumed triangular cross section curve, which was used the com- 
putation activation curve from the spectra Johns al. The points were recomputed for the 
cross section curve the photon difference method. 


IV. Recomputation Experimental Activation Curves 


The activation curves determined this laboratory for number photo- 
nuclear reactions and originally solved for the corresponding cross sections 
the total spectrum method have been recomputed the photon difference 
method. The activation curve for the reaction plotted 
Fig. This curve was originally measured Johns (13), but some addi- 
tional high energy points measured Mr. Baker have been included. 
The activation differences read from the above curve and the smooth curve 
fitted them are shown Fig. The calculated cross section curve plotted 
Fig. This and other computed cross section curves are believed 
reliable only the point where the cross section falls one-third the 
maximum value the high energy side the peak. the higher energy values 
were trusted, accuracy better would required for the entire 
activation curve. 


Because its widespread interest (see for example 17, 18, 27) the re- 
computed cross section curve for the reaction shown Fig. 
Haslam, Johns, and Horsley (8), using the total spectrum method, obtained 
very skewed cross section curve for this reaction. believe this due 
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is 16 18 20 22 24 «26 28 30 
MAXIMUM PHOTON (MEV.) 

Fic. activation curve for the reaction smooth curve has been drawn 
through the original points Johns al. plus some new higher energy points measured 
Baker. 


SATURATED ACTIVITY PER MOLE 
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Fic. activation difference curve for the reaction The points are the 
differences 0.5 Mev. intervals the smooth curve Fig. 
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SECTION 


CROSS 


PHOTON ENERGY (MEV) 


Fic. cross section curve for the reaction computed the photon 
difference m.thod, using the smoothed activation differences Fig. tail this curve 
very tnaccurately determined. 


PHOTON 


Fic. The cross section curve for the reaction recomputed the photon differ- 
ence method from the activation points Haslam, Johns, and Horsley. 
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the progressive smoothing characteristic the total spectrum method. The 
present calculation has yielded more symmetrical curve. Owing the narrow- 
ness the peak and the small energy difference between peak and threshold 
for the reaction, the computed tail for this curve probably somewhat more 
reliable than those for the other reactions reported here. 


The measurement the 2.55 minute activity resulting from the 
irradiation sulphur special interest because both the reactions 
d)P® and can contribute it. This activation curve has 
been measured more accurately Katz and Penfold (15) than that for any 
other reaction this laboratory. Any broad structural features the corres- 
ponding cross section curve would therefore expected show its solu- 
tion. The activation differences, shown Fig. 10, indicated slight bulge 
the lower energy side their peak, but the experimental points were 
accurately determined that this bulge could not smoothed out according 
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ACTIVATION DIFFERENCES MINUTE PER 0.5 MOLE 


19 20 2) 22 23 24 25 26 27 26 
MAXIMUM PHOTON ENERGY 


10. The activation difference curve for the reactions The points were 
taken from the original activation curve Katz and Penfold. The experimental accuracy estab- 
lishes the bulge the rising portion this curve. 
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the procedure recommended above. Therefore the first but not the second 
derivative the activation curve was smoothed this case. The resulting 
cross section curve plotted Fig. 11, with the low energy end magnified 
factors and 100. During this computation was discovered that the 
published cross section values (15) were too low about factor owing 
error the original computation. 


(MILLIBARNS) 


SECTION 


CROSS 


ENERGY (MEV) 


Fic. cross section curve for the reactions np)P® computed the photon 
difference method from the activation differences Fig. 10. The curve might represent super- 
position and np) peaked curves shown dotted the figure. 


The cross section curve Fig. might interpreted superposition 
two peaks, probably about indicated the dotted lines the 
figure. The lower these could identified with the reaction and 
the higher with the reaction The revised integrated cross section 
for the second these two reactions 0.013 Mev.-barns, about per cent 
integrated cross sections this region the periodic table. This 
the order magnitude predicted for this type reaction from neutron 
yields (5). 


also possible explain the two apparent peaks entirely due the 
(y, mp) reaction Following photon absorption, neutron evaporates from 
leaving various states excitation. The threshold for the emission 
proton from about Mev. lower than that for emission second 
neutron. Thus, the photon energy raised, the reaction sets in. 
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slightly higher photon energies the (y, cross section lowered owing 
competition with the (y, reaction, and rises again still higher photon 
energies where the protons can penetrate the barrier more efficiently. 
then not possible evaluate the contribution the cross section curve 


Fig. 
should noted that the bulge the activation difference curve did not 
cause oscillations the successively computed cross section values for the 


above reactions. 


Values the cross sections for the above four reactions are listed intervals 
Mev. y-ray energy Tables and III. Also listed are the recomputed 
cross sections for reactions Cu®, (13); (15); 
(14); and the reaction (9). each case the 
original experimental activation points served the starting point for the 
recomputations. The characteristic features the cross section curves 
originally determined the total spectrum method are listed Table IV. 
The redetermined quantities corresponding these computed the photon 
difference method are listed Table major changes were introduced 
the recomputation, but the newer values are believed more reliable. 
The new energies for the cross section maxima tend slightly higher than 
the old, but their differences near have been much reduced. The new 
cross section curves also appear somewhat smoother and more symmetrical 


than the old ones. 
The reaction was originally measured relative 


Cu®. Unfortunately the reaction contributes appreciable ac- 
tivity near the photonuclear threshold, and the lower half the cross section 


TABLE 


PHOTONUCLEAR CROSS SECTIONS IN MILLIBARNS FOR VARIOUS REACTIONS AT 1 MEv. 
INTERVALS AS REDETERMINED BY THE PHOTON DIFFERENCE METHOD 


Mev. 
14.4 104 150 
13.1 8.0 

8.8 6.0 

5.2 4.3 

3.3 

1.9 
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TABLE 


PHOTONUCLEAR CROSS SECTIONS IN MILLIBARNS FOR VARIOUS REACTIONS AT 1 MEV. 
INTERVALS AS REDETERMINED BY THE PHOTON DIFFERENCE METHOD 


energy 
Mev. 

130 

264 

470 

658 

665 

108 591 

386 

123 189 

116 

100 


123 
minimum 
140 
249 
348 
362 
313 
204 
100 
TABLE 


002 
0.015 
0.027 
0.040 
0.0185 
0.048 062 
0.071 
1.62 0.075 
2.35 0.075 
6.0 (0. 050) 


PREVIOUSLY PUBLISHED VALUES FOR THE CHARACTERISTICS OF VARIOUS PHOTONUCLEAR 


CROSS SECTION CURVES. 


THE SYMBOLS ARE USED THE TEXT 


Reaction 


n)Cu 


Cu®(y, 
n)Sb'” 


Em(Mev.) 


Maximum 


(barns) 


0.0116 
067 
060 
0.10 
>0. 
0.34 
>0. 078 


(0. 0004) 
(0.0015) 
00007 


Half-width 


SOA 


ans 


(Mev.-barns) (Mev.)*** 


Integrated 


3.7 
0.13 6.6 
4.5 
0.33 6.5 
6.6 
8.8 
0.83 6.9 
5.2 
2.0 5.2 
5.5 
(0. 0012) (4.9) 
4.4 
00041 


approximate lower peak the original cross section curve. 


**The principal cross section peak. 
difference between peak cross section and threshold. 
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q 
19.0 
18.3 
18.5 
17.5 
19.0 
18.5 
14.5 
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13.5 
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TABLE 


THE CHARACTERISTICS OF VARIOUS PHOTONUCLEAR CROSS SECTION CURVES AS REDETERMINED 
BY THE PHOTON DIFFERENCE METHOD. THE SYMBOLS ARE AS USED IN THE TEXT 


22.9 0.0131 2.8 0.046 4.2 
19.5 0.0167 6.5 129 
18.7 0.067 6.3 0.48 
18.5 0.054 5.6 0.34 6.5 
18.6 0.151 7.0 8.4 
18.7 124 7.9 0.99 
14.8 4.8 1.92 5.5 
4.5 


*Assuming the lower dotted curve Fig. 11. 
**Assuming the higher dotted curve Fig. 11. 
difference between peak cross section and threshold. 
the peak composed deuterons. 


curve undoubtedly distorted this effect. The correction the observed 
Cu® activities for self-absorption the samples was large, that the large 
absolute values the cross section are subject considerable error. Similar 
large values were obtained the University Pennsylvania (4). 


Some questions also arise connection with the reactions the 
antimony isotopes. single cross section curve was computed for the reactions 
both isotopes, since the ratios the activation points various X-ray 
energies were constant within experimental errors. Sb'” consists two isomers 
with half lives 3.5 minutes and 2.8 days. The shorter activity decays y-ray 
emission the longer, that measurement the 2.8 day activity gives the 
total yield formed irradiation. Its integrated cross section 1.92 
Mev.-barns agrees nicely with that predicted for antimony the basis 
photoneutron yields (5). Two isomers, min. and 6.0 day half lives, are 
reported for However, only the min. activity was measured. decays 
positron emission the end point the beta-spectrum being 1.70 
Mev. (3). When the yield this activity corrected for the recently published 
decay scheme (3), the integrated cross section for becomes 3.5 
Mev.-barns. This 80% higher than that predicted for antimony from photo- 
neutron yields (5). The decay scheme does not take into account small addi- 
tional positron component with end point 2.4 Mev., and the percentage 
K-capture competing with the 1.7 Mev. positrons assigned theoretical 
basis. The absolute yield therefore subject considerable error. 


One might expect, however, that the 6.0 day isomer should pro- 
duced yield comparable that the min. isomer. that case the inte- 
grated cross section for would very much more excess 


4 
ry 
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that predicted for antimony from photoneutron yields. The 6.0 day activity 
was observed Lindner and Perlman (20), who assigned from con- 
sideration complex such activities produced their irradiations Sn, 
Sb, and with high energy deuterons. decays K-capture and emits 1.1 
Mev. y-rays but beta-particles. therefore probable that the min. 
metastable state this isotope would about 1.7 Mev. above ground. 
very surprising that there should intervening states which the meta- 
stable level could decay y-ray emission. 


0.57 Mev. reported the literature accompany the decay 
the 2.8 day level During the course the antimony irradiations this 
laboratory lead absorption curve this activity was taken which indicated 
y-ray energy 0.58 Mev., good agreement with the above value. This 
experiment indicates that the 1.1 Mev. y-ray accompanied every disin- 
tegration the 6.0 day then the yield this isomer less than per 
cent the yield The authors therefore question the assignment the 
6.0 day activity and believe that the yield the min. activity 
represents the major portion not all the reaction 


When the cross section curve for (y, was recomputed from the 
original activation curve, large negative cross sections were obtained for photon 
energies above Mev. was therefore suspected that the activation curve 
was error. This curve has recently been remeasured Haslam, Smith, and 
Taylor (10), who found appreciably higher activities for betatron energies 
above Mev. The cross sections reported for this reaction Table III were 
therefore calculated from their new activation curve. The decay scheme for 
not known, and the values are therefore the cross sections, 
obtained correcting the activity observed the samples for geometry and 
self-absorption. They are undoubtedly much too low. 


Since the reaction was measured relative 
monitor, the (y, activation curve needed correction. This accounts for 
most the differences between the new and old cross section curves for this 
reaction. 


With the redetermined cross section characteristics listed Table 
possible re-examine the relation between the peak energy 
(Mev.) and mass number Steinwedel, Jensen, and Jensen (32) have pre- 
dicted variation whereas Goldhaber and Teller (7) prefer 
Fig. shows log-log plot the recomputed values E,, for 
reactions versus The best fit for the above type relationship given 


(14) 
This favors the prediction Goldhaber and Teller. However, the five middle 
points all lie above the straight line, whereas the five extreme points lie below 
Note added proof: Blaser, Boehm, Marmier, and (Helv. Phys. Acta, 24: 


1951.) report that they have been unable find the 6.0 day activity from irradiation 
antimony and proton bombardment tin. 


NED 
ge 
ar 
ns 
Ly 
rs 
re 
- 
O- 
O- 
e- 


536 CANADIAN JOURNAL PHYSICS. VOL. 


> 
z 
oO 

” 
” 


10 20 30 50 70 100 200 


21.8 


ENERGY MAXIMUM 


MASS NUMBER 
Fic. log-log plot the energy for the maximum cross section (Em) versus 


mass number. (b). semilog plot the same quantities. The equations the best straight lines 
are indicated the figure. 

it. This relationship therefore not too satisfactory. semilog plot the 
same quantities shown Fig. The best straight line given 


21.8 exp (—0.00268 A). (15) 


The constants this empirical equation are nearly the same those originally 
proposed (5) from computations the total spectrum method; the peak posi- 
tions have been raised about 0.3 Mev. This type relationship appears 
more satisfactory that the existing points distribute themselves about the 
straight line random manner. 


Summary 


new method has been proposed for the computation photonuclear cross 
section curves from their X-ray activation yield curves. based the 
photon differences between normalized Schiff bremsstrahlung curves. basic 
assumption made that measured activation points can fitted with curve 
having smooth first and second derivatives. empirical justification has been 
obtained for this smoothing procedure. The computational procedure simpler 
than that the total spectrum method Johns al., and most the arbi- 
trariness the solution removed. 


The method requires that samples irradiated behind cm. Lucite and 
that X-ray dosages measured Victoreen ionization chamber under the 
same conditions. However, little error introduced the irradiations are 
carried out under other conditions relative standard reaction whose acti- 
vation curve under the standard conditions known. also satisfactory 
obtain relative activation curve under other conditions, long doses are 
monitored Victoreen chamber Lucite, and normalize the curve 
standard irradiation. 
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Tables the Schiff curves P(E, Eo) and special functions and 
Eo) defined Equations (7) and (8) are appended this paper. 


Twelve activation curves measured this laboratory have been solved for 
the cross section curves the photon difference method. The results are 
contained Tables II, and The values differ little from those computed 
the total spectrum method, but the authors believe them more reliable. 


New constants have been determined for two suggested relationships between 
(y, peak energies and mass number, indicated Equations 
(14) and (15). 
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APPENDIX 


TABLE SCHIFF BREMSSTRAHLUNG SPECTRA MODIFIED CM. LUCITE ABSORBER AND 
NORMALIZED TOTAL INTEGRATED INTENSITY 100 ROENTGENS MEASURED 
VICTOREEN CHAMBER. VALUES GIVE PHOTONS PER CM.? PER 100 PER 1/2 MEV. INTERVAL. 
MULTIPLY FIGURES TABLE 


0.25 Mev. 
Photon 
energy, 

11.75 4.547 
12.25 3.486 
12.75 1.797 


APPENDIX I—Continued 


0.25 Mev. 


13.75 2.834) 4.285) 4.665) 4.766) 4.834 
2.656) 4.035) 4.249) 4.488 


538 
. 
Photon 
16.75 1.041) 1.897) 2.421 
17.25 1.776 
17.75 897 
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455 
10. 546 
8.793) 8.733 


7.445 
6.905 
6.417 
6.007 


4.913 
4.553) 4.598 
4.291 
3.984 


3.550) 3.672 
2.593 
2.139 


1.555 
0.780 


840 


APPENDIX I—Continued 


11.314 
10. 445 


0.25 Mev. 


19.75 20.25 20.75 


10. 343 


21.25 21.75 


22.25 22.75 


539 
Photon 
8.75 9.914) 9.801) 9.685 
13.75 4.952) 4.953) 4.944) 4.937 
15.75 3.757| 3.867 3.956) 3.965 
19.75 1.286) 2.246) 2.358 
21.75 1.341 
22.25 0.501; 0.967 
22.75 471 
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APPENDIX 
energy, 

16.75 3.567| 3.588) 3.594) 3.592) 3.582) 3.574 
20.25 2.247) 2.550) 2.689) 2.728) 2.766 
25.75 1.109) 1.264 
26.75 0.864 
27.25 0.604 
27.75 290 
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APPENDIX 


THIS TABLE MUST USED COMPUTE THE MODIFIED ACTIVATION DIFFERENCES BAa 


8.25 13.25 0.986 1.973 756 
8.75 0.551 13.75 1.052 2.125 23.75 3.975 
9.75 628 14.75 1.201 2.435 24.75 4.390 
10.75 0.710 15.75 1.379 20.75 25.75 4.820 
11.25 16. 21.25 2.927 26. 5.045 
11.75 0.813 16.75 1.591 21.75 26.75 260 
12. 865 17.25 1.710 22.25 3.307 27.25 5.495 
12.75 17.75 22.75 3.516 27.75 5.720 


) 
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APPENDIX III 


TABLE Eo). THE SUCCESSIVELY DETERMINED CROSS SECTIONS 
MUST MULTIPLIED ALONG THE HORIZONTAL ROWS THIS TABLE DETERMINE THE 
QUANTITIES SUBTRACTED FROM THE MODIFIED ACTIVATION DIFFERENCES 


Photon 0.25 Mev. 
energy, 

8.25 0.988 0.885 0.621 0.398 0.268 0.172 0.135 0.090 0.046 
8.75 0.967 0.802 0.553 0.362 0.244 0.170 0.130 0.077 
9.25 0.949 0.713 0.485 0.325 0.222 0.167 0.122 
9.75 0.935 0.630 0.440 0.294 0.213 0.162 

10. 0.929 0.562 0.419 0.288 0.210 
10.75 0.916 0.486 0.394 0.286 
0.903 0.465 0.386 
11.75 0.882 0.455 
872 


APPENDIX 


Photon 0.25 Mev. 
energy, 


14.75 0.856} 0.493) 0.149 
0.584 


| 
| 
| 
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APPENDIX 


Photon 0.25 Mev. 


18.25 18.75 19.25 19.75 20.25 20.75 21.25 21.75 22.25 22.75 


0.068 


19.75 0.950} 0.441) 0.344) 0.238 
21.25 0.465 

0.678 


E | | | | | | 
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544 
APPENDIX 
Photon Eo 0.25 Mev. 
energy, 
23.25 23.75 24.25 24.75 25.25 25.75 26.25 26.75 27.25 27.75 

12.75 —0. —0. —0. —0. —0. —0. —0. —0. —0. 096 
13.75 —0. —0. —0. —0. —0. —0. —0. —0. —0. 072 
0.766) 0.514) 0.465 
25.75 0.593) 0.535 
26.75 0.810 
27.25 1.043 
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Examples electron events. 


(a) Electron pair low quantum energy. 


(b) The 


ation energetic photon close and parallel with the track 


particle. reasonable suppose that the photon represents Bremsstrahlung 


from the parent 
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ELECTROMAGNETIC CASCADES PHOTOGRAPHIC 


Abstract 


Photographic emulsion techniques are described for observation and identi- 
fication electromagnetic processes which occur the soft component the 
cosmic radiation. Several cascade showers are illustrated and measurements 
these are discussed the light theoretical predictions. 


Introduction 


The introduction photographic emulsions sensitive the passage 
relativistic particles charge /e/ has led new results the study electro- 
magnetic processes initiated the cosmic radiation. number experiments 
designed investigate the properties the fluxes energetic electrons and 
photons great altitudes have been carried out this laboratory during the 
past two years. Ilford emulsions 400 thick and capable recording the 
tracks electrons ionizing the minimum rate were carried great altitudes 
free balloons filled with hydrogen. After recovery and suitable photographic 
processing, the emulsions were examined systematically through microscopes 
fitted with oil immersion objectives. 


the course this search, sets nearly parallel tracks low and constant 
grain density have been observed. When the points origin these tracks 
can found the emulsion, usually possible interpret them electro- 
magnetic events one the following types: 


(i) Electron* pair production y-ray. (Fig. 
(ii) Direct pair production fast electron (‘trident’). 
(iii) The apparent conversion Bremsstrahlung y-ray electron pair 
beside the track fast electron (Fig. 10). 
(iv) Energetic electron recoil from fast particle. 
(v) Multiple electron pair production which four electrons are apparently 
created the conversion single y-ray (12). 
(vi) Cascade generation which may involve any all the above processes. 
Because the emulsion thin, cascade events complete enough provide 
material for analysis are rare. The observation and measurement these 
phenomena the emulsion permits insight into cascade development 
‘hitherto unattainable. 


Manuscript received May 30, 1951. 
Contribution from the Wills Physical Laboratory, University Bristol, Bristol, 


England. 
The term ‘electron’ used here designate both positrons and negatrons. 
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Techniques Observation 


Experience has shown that the X45 objective suitable for the arduous and 
prolonged search which involved the accumulation large numbers 
these events. has been necessary limit the observations those events 
path length the emulsion sufficient permit energy determination 
multiple scattering (9). 


The following methods search have therefore been adopted. Parallel strips, 
across the plate, 250 wide, are examined closely. These are spaced dis- 
tance that will ensure the discovery the majority events, path length 
the emulsion exceeding the required value. This procedure greatly increases 
the rate which significant events are detected. 


The majority such events found are single electron pairs produced 
y-radiation. pair most readily detected two apparently related and 
nearly parallel tracks minimum grain density, and rarely found through 
the direct perception its point origin. systematic search marked 
discrimination may noted against finding those pairs for which the paths 
the two particles are not appreciably parallel for any distance the emul- 
sion. This reduction search efficiency applies those pairs low quantum 
energy, those marked disparity the energies the two electrons, and 
those energetic pairs which the two electron tracks never become notice- 
ably separated the emulsion. 


Further difficulties discrimination arise with regard the selection 
material for measurement. has been found convenient apply criterion 
minimum path length the emulsion for each track event, but this 
not necessarily the most desirable for analysis data obtained multiple 
scattering. Thus, few pairs low quantum energy satisfy length require- 
ment more than mm. per track, whereas the errors measurement 
the tracks very energetic pair particles may disproportionately high for 
path lengths limited only mm. the emulsion. The length criterion may 
also cause the rejection pair marked disparity the loss the more 
scattered track through one the emulsion boundaries. The same effect ap- 
plies those pairs wide angle which separate plane nearly perpendicular 
that the plate. 


Another source experimental error arises these observations from varia- 
tions the efficiency individual observers. Such factors the conditions 
physical and mental comfort prevailing the laboratory may noticeably 
influence the data. Visual faults may also exercise some effect, such possible 
discrimination astigmatic observer against pairs orientated par- 
ticular angle the field view. These factors vary from day day and from 
observer observer and are always difficult assess. 


Identification Electron Tracts 


The track electron with total energy exceeding 1.5 Mev. exhibits 
specific ionization within 10% the minimum value. Since the minimum 


q 
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specific ionization relativistic particles charge /e/ independent mass, 
further identification the tracks energetic electrons must made 
one more the following characteristics. 


(i) Multiple Scattering 


mean angle scattering greater than 0.25° per 100 path length indicates 
momentum less than 100 Mev./c for the particle. This identifies the par- 
ticle electron since the produce tracks minimum ionization 
only when the mean angle scattering less than 0.25° per 100 as- 
sumed that charged particles exist intermediate mass between electrons 
and 


(ii) Radiative Energy Loss 

The theory Heitler shows that the cross section for the energy loss 
radiative collisions proportional the reciprocal the square the mass 
the incident particle. Therefore the probability very small that the process 
will found occur the tracks particles other than electrons. sudden 
increase the mean angle scattering track thus indication the 
electronic nature the particle. Apparent energy loss radiation has fre- 
quently been detected the tracks particles originating pair events but 
has never been found occur the tracks particles emerging from nuclear 
explosions. The work Camerini (7) has established that electrons are 
rarely, ever, emitted the course nuclear disintegrations. 


(iii) Phenomenological Evidence 


The foregoing criteria provide evidence that the majority particles origin- 
ating such events tridents, multiple pair production, and cascade pro- 
cesses, are electronic mass. addition reasonable assume that par- 
ticles more massive than electrons are not created the materialization 
y-rays. When identification based these phenomenological grounds (Fig. 
1b), must remembered that cascade showers may possibly form close to, 
and direction nearly parallel with, the track relativistic meson 
proton. This situation arises through the simultaneous generation charged 
and neutral energetic nuclear explosions. The majority cascade 
showers great altitudes originate from the conversion the photons which 
result from the spontaneous decay the short-lived neutral 


Characteristics Electron Tracks 
(i) Loss 


Bloch (5) first derived the relation between the rate energy loss 
ionization, and the velocity, for charged particles. calculation, assuming 
perfect screening the nucleus the atomic electrons and mean value 
the ionization potential 13.5 ev. per atom, yields 0.57 Mev. the minimum 
ionization loss per mm. charged particle the NT4 emulsions. 
The order magnitude this result has been confirmed the experiments 
Brown al. (6). 
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(ii) Radiation Length 

The radiation length, Xo, defined the path length which the initial 
energy fast electron reduced through radiative collisions, the average 
factor (3). For great energies, which this rate energy loss nearly 
constant, given the equation:— 


where the number atoms per cc. the medium and the 


the electron. Evaluation for the nuclear 


classical radius 


emulsion yields value 2.9 cm. Measurements have been made the mean 
angle scattering, over mm. intervals electron tracks exceeding 1.5 cm. 
the emulsion. Preliminary results shown Fig. verify the order magni- 
tude the theoretical result. The ordinates indicate the mean value the 
ratio each interval distance along the tracks. interest 
note that the method independent the scattering constant and the 
initial energy the particles. 


X in mm. 


Fic. radiative energy loss tracks fast electrons. The slope the 
curve yields the value 3.4 cm. for the radiation length the emulsion. The limits error indicated 
are the standard deviations. 


Critical Energy 

When fast electrons traverse matter, energy lost both ionization and 
radiative collisions. The rate loss the latter process varies with energy, 
and the particular energy which the radiative loss equal the ionization 
loss called the critical energy. may noted that radiation length, 
electron loses, through ionization, energy nearly equal the critical energy. 
This quantity characteristic the absorbing material, and for nuclear emul- 
sions the neighborhood Mev. 


(iv) Conversion Length 
The conversion length that which beam photons reduced pair 
production the initial intensity. This length decreases with increasing 
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quantum energy nearly constant value above 300 Mev. From theory 
may seen that the conversion length the emulsion great energies 
3.8 cm. 9/7 the radiation length. 


Cascade Showers 


the course the examination about cc. electron-sensitive photo- 
graphic emulsion exposed the cosmic radiation 68,000 ft., eight striking 
examples electron cascade phenomena have been found. Facsimile drawings 
representing the tracks three these multiplicative events are shown 
Figs. and similar example has been published previously (13). 
these drawings the lateral scale has been expanded factor greater 
relative the axial scale, since otherwise the high collimation the electron 
tracks would obscure the details the cascade development. 


The cascade process was first investigated theoretically Bhabha and 
Heitler (4) and Carlson and Oppenheimer (8). Experimental measurements 
Wilson cloud chambers have been made many workers. Typical these 
investigations are those who studied the development showers 
absorbers various thicknesses. Others, such Lovati al. (15) have 
applied cascade theory the determination great quantum energies the 
course investigations the 


The photographic method offers peculiar advantages the study these 
‘soft’ showers. The points origin and the spatial distribution the indi- 
vidual events are readily determined relative the cascade axis. This enables 
reasonable discrimination between tridents and electron pairs arising from 
the conversion Bremsstrahlung. Wilson cloud chambers, the details 
the multiplicative processes are necessity obscured within the thickness 
metal absorbers, whereas these studies the sensitive material itself the 
absorber. The degeneration the initial energy among the daughter particles 
may analyzed all stages the methods multiple scattering. 


These events are rare because the primary must high energy and 
closely parallel with the plane thin emulsion, that the secondary tracks 
have reasonable probability remaining within the sensitive layer. cor- 
rection has made the observations take into account continuous 
escape low energy particles through the emulsion boundaries even the 
initial stages the cascade. anticipated that will prove possible within 
the next few years manufacture, process, and examine considerably thicker 
emulsions. Such improvements will greatly minimize these difficulties. 


our study cascade showers have compared our observations with 
deductions from the one-dimensional theory tabulated Janossy (14). 
Following this author have employed the quantity defined Equation 
the unit cascade length. should noted that this not the same defini- 
tion cascade length that used some other workers, for example 
Arley (1). Since the efficiency detecting those electrons less than Mev. 
cascade poor, have restricted our computations the numbers 
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10 


Fic. facsimile drawing representing the tracks cascade shower observed the emul- 
sion. The shower originates the conversion energetic photon. The plain and dot-dash 
lines indicate electrons pairs produced The dashed lines indicate the tracks two 
secondary electrons low energy from the trident. Lateral scale 20. 
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Fic. reproduction the tracks soft shower originated fast electron (plain line). 
The dashed lines indicate trident secondaries and the dot-dash lines the tracks electrons pairs 
which are produced Bremsstrahlung. Lateral scale 20. 
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Fic. electron (full line track) initiates cascade one trident (dashed lines) and ten 
electron pairs (dot-dash lines). Lateral scale 50. 
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electrons expected above this energy. Also, because the loss energy above 
Mev., the critical energy, predominantly radiation, have neglected 
the ionization term the cascade diffusion equations. can then shown 
that the number electrons photons, with energies exceeding pro- 
duced primary energy after length cascade units given the 
integral: 


where the particular functions and depend whether the cascade was 
initiated photon electron. Further details may found Janossy’s 
book (14). This author also gives tables derived the use the saddle- 
point method approximating the integral. Since the showers which have 
found the photographic emulsions are confined observable lengths the 
order only one two cascade. units, necessary retain the second 
exponential Equation 


Observations the Illustrated Cascades 


example the application Equation the observed showers, 
Fig. discussed some detail. This drawing represents cascade initiated 
addition, associated electron pairs are formed the conversion 
strahlung photons. All the tracks ultimately disappear result the elec- 
trons passing through the air glass interfaces the emulsion. The origin 
the first pair taken the beginning the shower, because the exact 
distance travelled the primary photon before conversion not known. The 
shower now considered the sum two individual cascades, the primaries 
which are the two electrons the first pair. The energies these are 10,000 
and 8000 Mev. this basis, the average number electrons energy above 
Mev. expected for 4.7 and 4.4 respectively, total 9.1. For 
corresponding numbers are 13.9 and 12.4, which give about 26. The 
initial energies all electrons observed the shower are set out Table 
The experimental number any value taken made two 


parts: 
(a) those tracks actually seen; and 
(b) those electrons which were created earlier, left the emulsion before 
and which probably would have had energy exceeding Mev. 
Taking account these considerations, find for the shower under consider- 
ation that for values and the observed numbers are and 
respectively. 
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TABLE 


Event Mev. Mev. 
Pair 10,000 8000 
Pair 2700 1400 
Trident (7900 Mev.)* 
Pair 475 
Pair 
Pair 
Pair 425 300 
Pair 180 
Pair 130 
Pair 135 
Pair 260 
Pair 


trident has three secondary particles. 


should pointed out that cascade theory predicts expected numbers 
only, and that fluctuations the observed values from the mean may very 
large (16). Another factor which influences the observations the failure 
observe cascade development which would occur the emulsion were not 
limited thickness. This effect difficult estimate but depends upon the 
initial energy and the value considered. The relatively small number 
factors. 


The majority the cascade showers found far arise, however, from the 
degradation energy one more fast electrons which enter the emulsion. 
conceivable that under these circumstances one more photons may 
accompany the electrons. least one shower clear that this has fact 
occurred since secondary pair has been formed prior the entry from the 
glass the parent electron. 


Another factor which appears importance the development 
cascades the trident process (13). This modifies the rate production 
secondary electrons great energies and effect which has, the 
present, been neglected the diffusion equations. have found that five out 
seven the showers incorporate one more tridents. should noted, 
however, that account unequal energy division among the secondary 
particles, this process relatively less important than would equi- 
partition energy invariably occurred. From our earlier studies have 
observed that, more than 90% the tridents, more than half the available 
energy given one the secondary electrons. Indeed, rather more than 
50% the examples, one the outgoing electrons carries more than two-thirds 
the primary energy. 


The cascade reproduced Fig. interesting example soft shower. 
this event three tridents occur, together with four associated pairs, 
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path length 1.5 cascade units. The energy the primary electron about 
1.3 Mev.; this energy, the cross section for trident production ap- 
proximately 1/10th that for radiative collisions (17). 


this cascade, the development begins with trident, and this followed 
the creation pairs path Jength 4.5 cm. The number electrons 
mental number estimated 12. this event the geometrical loss probably 
small, because the event remains close the median plane between surface 
and glass. The large discrepancy observed may explicable terms any 
the other effects mentioned above. 


Origins Cascade Primaries 


Considerations the decay processes charged and neutral 
indicate that altitudes exceeding km. the great majority electron 
cascades originate from the y-rays which occur through the decay the short 
lived particles. Some these y-rays may materialize the emulsion, re- 
sulting energetic cascade (Fig. 3), while others may materialize the 
atmosphere considerable distance from the plate assembly, that the 
observations pair associated electrons entering the plates would rare. 
The single electrons arising such circumstances are considered form the 
majority our shower primaries. (Figs. and and others not illustrated.) 
Rossi (18, 19) has shown that only very small proportion fast electrons 
entering the emulsion these altitudes can due the spontaneous decay 
energetic recoil processes. the emulsion, however, the 
probability energetic electron recoils from relativistic mesons much greater. 
Two cascades arising this way have been identified this laboratory; one 
these has been previously published (13). 


Conclusions 


The observations which have been described lead the conclusion that 
the photographic method provides new and powerful approach the de- 
tailed study electromagnetic shower phenomena. expected that present 
limitations will largely removed when becomes possible process thicker 
emulsions without distortion and examine these with reflecting micro- 


scopes (2). 


predicted the mathematical theory that there are very large statis- 
tical fluctuations (16) the numbers particles expected any given stage 
‘in the development cascade. the light this result, the present obser- 
vations are not inconsistent with the cascade theory. have also been able 
examine detail the individual processes cascade formation. par- 
ticular, the importance the trident process the initial stages energetic 
cascades has been demonstrated. 
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200 HIGH TENSION SET FOR THE ACCELERATION 
AND He* 


Abstract 


200 kv. high tension set suitable for the acceleration and de- 
scribed. Special features incorporated the accelerator include 
frequency type ion source capable giving large beam currents with low gas 
consumption, differential pumping system for recovering the gas used the 
ion source, and electronic stabilizer which reduces voltage fluctuations 
negligible value. 


Introduction 


The advantages using and bombarding particles the study 
nuclear reactions have Jong been recognized. Both nuclei have low binding 
energies per nucleon and consequently should induce reactions that are 
general exothermic and that, for light elements, might expected have 
adequate yields low energies. addition compound nuclei are formed which 
have many possible modes disintegration some cases these lead product 
nuclei that cannot easily studied any other means. 


Until recently neither nor has been available sufficient quantities 
allow systematic study their interactions with other nuclei. 
emitter, decaying with half-life 12.5 yr. into which, although 
stable, comprises only atmospheric helium. However, the advent 
high power nuclear reactors has made possible, through the neutron bombard- 
ment lithium, the production sufficient quantities tritium and hence 
begin such systematic study. The purpose this paper describe 
200 kv. accelerator which has been built Chalk River for the bombardment 
light elements with tritons and ions. Three special features are incor- 
porated the accelerator, the first two being dictated the scarcity gas 
available and the radioactivity These are: 

(1) radio-frequency type ion source capable producing beam currents 
the order 200 with hydrogen gas consumption cc. per hr. 

(2) differential pumping system with mercury diffusion pumps which 
allows recovery more than 95% the gas admitted the ion source; 

(3) electronic stabilizer which reduces voltage fluctuations Jess than 
0.05%. 

General Description 


The main features the high tension set are shown Fig. The accel- 
erating potential provided two transformer-rectifier units, and 


Contribution from the Nuclear Physics Branch, Atomic Energy Project, National Research 
Council Canada, Chalk River, Ont., Canada. Issued N.R.C. No. 
Physicist, United Kingdom Staff; now the University Liverpool, England. 
Now with Schlomberger Well Surveying Corp., Houston, Texas. 
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the front the accelerator. 


General view 
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cascade. Power fed these through the isolating transformer The large 
structure the power stack. metal base covered with plywood supports 
insulating column, ft. long and ft. diameter made Grade Rolled 
Dilecto Tubing chosen for its combined strength and high resistivity. top 
this the high tension enclosure, wood and masonite box having smooth 
exterior painted with conducting paint* provide equipotential surface. 
This houses kva., 500 cycle generator providing power for diffusion pumps 
and all electronic equipment high potential. The generator driven 
Jong belts from h.p. repulsion induction motor the base and carbon- 
pile regulated. Extending horizontally from one wall the high tension box 
cylindrical aluminum tunnel ft. long and ft. diameter, capped 
the far end ft. diameter hemispherical spinning (Fig. 2). This tunnel 
contains the ion source, differential pumping system, and gas recovery reser- 
voirs. The upper half the made removable wire mesh allow 
observation and access the ion source. Directly below the ion source may 
seen the accelerating conventional design with two accelerating 
gaps. Two in. diameter Pyrex cylinders each ft. long form the vacuum 
housing. the bottom the tube the pumping manifold liquid air trap 
water baffle and diffusion pumps and backing pump O-ring vacuum 
valve (3) may used isolate the target assembly below. many 
experiments the beam, after passing through the Sylphon bellows bent 
through 90° resolving magnet into the target chamber. The H.T. voltage 
measured potentiometer method the grounded end the 250 megohm 
precision resistor chain part which shown Fig. 


Detailed Description 
(a) 200 kv. Power Supply 


circuit diagram the unit which produces the main accelerating potential 
shown Fig. consists two cascaded half-wave voltage doublers 
and Fig. 1), each capable producing ma. 100 kv. The chassis 
the two units are respectively kv. and 150 kv. above ground that 
necessary feed power them through isolating transformer The 
voltage varied means Variac the primary the isolating trans- 
former. series with the primary are various time delays and safety switches 
protect the high voltage supply. The ripple approximately 0.3% per ma. 
and reduced considerably the stabilizing circuit described below. Con- 
nections the accelerator are made the 100 kv. and 200 kv. levels through 
4-megohm 400 protective resistors. 


(b) Voltage Measurements and Stabilization 


The accelerating voltage measured with 250 megohm resistor chain 
(Fig. consisting XR-5 wirewound unitst megohms each with 
tolerance 0.1%. wa. meter, accurate 0.5%, arranged dividing 


Rubalt supplied Alfred Hague and Co. Inc., St., Brooklyn, N.Y., U.S.A. 
Supplied Shallcross Mfg. Co., Collingdale, Pa. 
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circuit the bottom the chain read full scale 300 kv. Provision made 
for inserting potentiometer place the meter for more accurate readings. 


schematic diagram the method stabilization shown Fig. 
has been fully described elsewhere (2). Briefly, type 4E27 stabilizing tube, 
capable handling variations several thousand volts, placed series 
with the high tension supply and its plate-to-cathode voltage controlled 
amplified error signal taken from the upper end the resistance chain. 
stabilization factor about achieved, which suffices reduce output 
voltage variations less than 0.05%. This, conjunction with constant 
current supply for the beam deflecting magnet, permits attainment very 
steady resolved beams. 


4E27 STABILIZER 
TUBE, 


O- 200 KV (STABILIZED) 


FEEDBACK 4 
AMPLIFIER 


4 MEGOHM 
SURGE RESISTORS 


250 MEGOHM 
HIGH STABILITY 
BLEEDER 


Fic. method stabilization the H.T. voltage. 


evident from Fig. that the resistor chain may connected either 
the end the focusing supply. The former connection has 
the advantages that the resistor chain measures the over-all voltage (with the 
exception 2-3 kv. extracting potential) and that this voltage stabilized. 
additional kv. may obtained using the second connection the 
expense some voltage stability. This done the study reactions with 
low yield order obtain maximum intensity emitted particles. 


(c) Ion Source 


Details the ion source are shown Fig. the radio-frequency 
type described Thonemann al. (5) and was chdsen for its ability pro- 
duce large beam currents high atomic percentage with low gas consumption 
and with small energy spread the emerging ions. 


this source Mc. per second r-f. discharge maintained air 
cooled pyrex envelope, gas being allowed enter the top. The flow hydro- 
gen tritium regulated palladium Jeak; when helium accelerated 
this leak replaced fine control needle valve.* order maintain 
high percentage atomic ions the glass envelope must cleaned occasionally 
and hence made easily removable. Hydrochloric acid very effective 
removing metal sputtered the inside walls. 


Obtained from Edwards and Co., Worsley Bridge Road, Lower Sydenham, London, S.E. 
26, England. 
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Fic. Details the r-f. ion source assembly. 


Ions are drawn from the discharge small dural probe containing 
canal 0.070 in. diameter and 3/4 in. long. This probe also made easily 
replaceable becomes charred and pitted after several weeks operation. 
was found necessary provide baffle plate below the probe prevent 
return electrons welding the support. 


The geometry the probe and the glass chimney surrounding fairly 
critical. Poor ion extraction results the height the glass chimney above 
the probe too great. the chimney too small diameter will suffer 
severe overheating from ion bombardment, while too large excessive 
probe current will drawn. The dimensions shown Fig. represent satis- 
factory compromise. 
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potential 0-5 kv. applied between the probe and top the ion source 
extracts ions through the probe, and further potential 0-40 kv. focuses 
the ions through second small canal, 1/4 in. diameter and in. long, 
into the accelerating tube. The lower electrode the focusing lens rests 
three cams attached shafts which protrude through the walls; these permit 
tilting the electrode secure maximum beam current through the second hole. 
the assembly the electrodes the two canals are aligned along the axis 
the accelerating tube. Design the r-f. oscillator and voltage supplies (Fig. 
quite conventional. However, the r-f. coil wound about the glass envelope 
must the same potential the extracting probe sparking through the 
glass will occur. The position the coil should adjusted give maximum 
beam current. not very critical. 


Typical operation the ion source follows. With 3.5 kv. extracting 
voltage, ma. extracting current, and kv. focusing voltage, target cur- 
rent 200 obtained with hydrogen gas consumption cc. per hr. 
N.T.P. This current the form unresolved beam collimated 1/8 
in. diameter. The beam contains about 70% monatomic, 10% diatomic, 
and 20% triatomic hydrogen. 


(d) Differential Pumping and Gas Recovery System 

More than 95% the gas admitted the ion source recovered dif- 
ferential pumping system (Fig. which maintains the pressure the focusing 
region about mm. mercury. This all metal system consists 
three mercury diffusion pumps, D.P.I. MHG-50 backed Edwards 
Type followed another pump the same type which has been modified 
operating against higher pressures, mercury pumps are preferred oil pumps 
because greater cleanliness and because hot oil may cause tritium 
exchange. The pumps are cooled transformer oil circulated via gear 
pumps and long Saran tubes through water cooled heat exchanger ground 


potential. 


Gas passing through the pumps accumulates cylindrical steel tank having 
large piston fitted with O-ring make vacuum seal with the wall the 
tank. Thence the gas forced into storage reservoir containing spongy 
uranium which absorbs tritium. the end run all the tritium absorbed 
the spongy uranium and the recovery system isolated from the main accel- 
erating tube means vacuum valve (3) prevent mercury contam- 
ination the ion source. the beginning new run the whole recovery 
system first pumped out with the mechanical backing pump. Heating the 
uranium then forces the hydrogen-tritium mixture into the 200 cc. bottle 
beside the ion source. This flask holds sufficient gas for many hours oper- 


ation. 


Spongy uranium excellent purifier well storage reservoir for tri- 
tium absorbs hydrogen isotopes and impurities low temperatures and 


. 


AlddNS Y3IMOd ONISNIOS AX 


VOL. 


00S 


HHH MOS 


3Y1N39 


MIOIAOX b 


AlddNS 
Y3IMOd ONILIVYLXI ANS-O 


CANADIAN JOURNAL PHYSICS. 


564 
z 
=On 
wr 
= >2 
HHH 
= 
= 
Qt 
N 
> 
° 
o 
Oa 


Q 
=x 
& 
=x 
NX 
w 
x 


ssv19 


SSV19 


SdWNd 


39NVS 
“ww O92 -O 


YIOAN3S34Y 


39NVS 
“ww 


WOINVHOAW OL 


\ 
Per 
565 
qo 
9 
oO 
\ 
w 
ocx 
or 
| 
oa 
ad 
oc 
Or 


566 CANADIAN JOURNAL PHYSICS. 


gives only the hydrogen isotopes when heated 450°C. comprehensive 
account the purification and manipulation tritium this system 
published elsewhere. 


the acceleration the recovery reservoir merely empty 1500 cc. 
container into which the gas forced the piston arrangement. The dead 
space the connection the piston assembly about one cubic centimeter. 


(e) Accelerating Tube and Main Pumping System 

The main accelerating electrodes are three inches diameter and are made 
highly polished mild steel, the two accelerating gaps being approximately 
1-1/2 in. length. These electrodes are centered screw adjust- 
ments using O-ring seals. Experience has shown disadvantage this type 
accelerating gap collection charge the glass and resultant shifting 
the ion beam. Re-entrant electrodes would preferable. Normally this 
shifting not serious and compensated applying variable voltage across 
pair parallel plates installed inside the vacuum manifold. This provides 
adequate compensation one plane while adjustment the magnetic field 
the resolving magnet provides compensation the plane right angles. Thus 
complete control the ion beam attainable. 


The main pumping system MC-275 oil diffusion pump 
using Octoil, equipped with both water baffle and liquid air trap. This 
D.P.I. MB-100 vertical booster using butyl sebecate. The fore 
pump Kinney Type 556. Pressure the pumping manifold measured 
with ionization gauge, operating pressures the order mm. 
mercury being normally obtained with beam on. arrangement Pirani 
gauge with Sensitrol meter automatically switches off the diffusion pumps 
the pressure the accelerating tube rises above 


(f) Controls and Metering 

Control the accelerator conforms standard practice and requires little 
elaboration. the power supplies and variable gas leak for the ion source 
and the high tension voltage stabilizer are high potential they must 
remotely controlled: Inside the Dilecto column are long bakelite shafts rotated 
selsyns the bottom, which turn are activated mate selsyns the 
control panel. switching effected means solenoid relays the 
base the power stack connected via nonconducting cord microswitches 
the top. 


panel the high tension stack with Lucite face (Fig. houses meters 
displaying the generator output voltage, r-f. oscillator current, extracting 
current and voltage, focusing potentia!, voltage across the stabilizing tube, 
and pressure the differential pumping manifold. Philips gauge used 
indicate this pressure. Two dial gauges and thermocouple gauge meter, which 
monitor various pressures the recovery system, are mounted the 
visible from the floor. 
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(g) Tritium Monitoring* 


safeguard personnel has been necessary provide means detecting 
accidental escape tritium into the laboratory. this purpose large exhaust 
fan and monitoring system with alarm bell have been installed the roof 
above the accelerator. This system consists two identical air ion chambers, 
one sealed off and one open the ends receive the air which proceeds from 
the laboratory via ion collector. The outputs the two ion chambers are 
connected d-c. amplifier that zero signal results the ion current 
equal the two chambers thus cancelling background due penetrating 
radiations. Ionization produced tritium decay the chamber through 
which laboratory air passing will produce signal which amplified ring 
alarm bell and open louvres thus increasing factor the rate 
exhaust air from the laboratory. 


Operation 


The high tension set described this paper has been operating for approxi- 
mately one year and has given very little trouble. This may attributed 
the effectiveness O-rings vacuum seals and the robustness the elec- 
trical supplies both for the main accelerating potential and the ion source. 


far the accelerator has been used chiefly for the study nuclear reactions 
induced tritons bombarding tritium and lithium targets (1, 4). the pre- 
sent time acceleration pure being undertaken and, He’ scarce, 
provision has been made for recovering the small amount gas passing down 
the accelerating tube. Swedish General Electric Type JHF33 mercury dif- 
fusion pump with rated backing pressure mm. mercury has been 
connected between the main diffusion pump and the Kinney backing pump. 
The set can operated for many hours with the Kinney pump isolated and 
all gas from the main tube being pumped through the oil and mercury 
diffusion pumps into liter reservoir. the end the run, the gas ac- 
cumulated the reservoir pumped back the accelerating tube the re- 
covery system. 


Whenever tritium accelerated a-particles from the very prolific reaction 


appear, even though deuterium gas has never been allowed the system. 
The natural isotopic content deuterium present hydrogenous oil vapor 
deposited the target sufficient give large yields a-particles pre- 
cautions are not taken. These a-particles are somewhat useful for calibration 
purposes but other cases constitute very annoying background. the 
acceleration was necessary change the ion source completely 
eliminate tritium contamination. 
Designed Carmichael and Hitchcock this laboratory. 
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SOME STUDIES ANGULAR CORRELATION! 


Abstract 


The beta-gamma-angular correlation has been measured function 
beta energy using channel kicksorter and thick crystal counter beta- 
detector. The differential correlation has been found change 
from Mev. —0.44 the end the beta spectrum. When inte- 
grated numerically over all beta energies greater than Mev. the integrated 
angular correlation coefficient obtained agrees with the value measured 
directly. Experiments the gamma-gamma-angular correlation and 
performed with the same apparatus are agreement with previous results 
other workers. 


Introduction 


The main purpose the present paper describe some experiments 
the angular correlation between the beta particles and gamma photons emitted 
cascade the Sb™ nucleus. particular the dependence the angular 
correlation the energy the beta particles was the object this detailed 
investigation, the results which have already been discussed briefly 
preliminary report (3).* 

The present experiments beta-gamma-angular correlation were 
undertaken because the results the subject available the literature were 
contradictory. Ridgway (10) reported the correlation coefficient asa —0.17, 
while the results Beyster and Wiedenbeck (1) gave —0.26. was 
rightly hoped that the integrated angular correlation was strong re- 
ported Beyster and Wiedenbeck (1), measurement the dependence 
upon beta energy could attempted, such energy dependence being theo- 
retically predicted Falkoff and Uhlenbeck (4), but far case being 
observed experimentally quantitative way. 


Prior the experiments beta-gamma-angular correlation some experi- 
ments the gamma-gamma-angular correlation and were done. 
Although the results these latter experiments are not new, Brady and 
Deutsch (2) having observed this previously, seems worthwhile discuss 
them briefly well. This done because the same apparatus was employed 
the beta-gamma experiments, the consistency the gamma-gamma 
correlation results lends support the subsequent beta-gamma-correlation 
results. Moreover, the accuracy obtained the present gamma-gamma cor- 
relation experiments perhaps slightly better than the previously reported 
work. 


received July 20, 1951. 
Contribution from the Department Physics, the University British Columbia, Van- 
couver, B.C. 

Although many nuclei have been investigated for correlation, only four 
cases, namely and has any correlation been observed. The results different 
observers for are contradictory, while confirmation the results for has been 
reported. The only report attempted measurement the beta-gamma-angular correlation 
function the beta energy very brief one Stevenson (11), which made use 
lens spectrometer. 
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Consequently, Part and Part the present paper will devoted the 
gamma-gamma- and the beta-gamma-angular correlation 
spectively. 


Gamma-Gamma-Angular Correlation Co® and 


The apparatus employed measuring the gamma-gamma coincidence rates 
function employed RCA-5819 photomultipliers conjunction with 
anthracene crystals gamma counters. The source was mounted 
the center spectrometer table with the counters placed about in. from 
it. The source was the form metallic cobalt, while the source was 
the form scandium oxide. Pulses from the photomultiplier were amplified 
single stage head amplifier and transmitted cathode follower 
linear pulse amplifier with rise time 0.15 and gain 100. 
Following the linear amplifier was pulse height discriminator and coin- 
cidence mixer with resolving time 0.17 The discriminator had dead time 
counting rates 5000 per sec. could handled with counting 
loss about 1%. With source strength selected give this counting rate, 
the accidental coincidence rate was about one half the true coincidence rate. 


The main difficulty experienced the coincidence measurements due 
radiation being scattered from one counter the other. The gamma rays 
interact with the counter crystal mainly Compton scattering. The scattered 
photon may enter the second counter and counted, thus producing false 
coincidence count. the counters have about the same efficiency for detection 
the scattered radiation for the emitted gammas, then these scattered 
coincidence counts may amount much 50% the true coincidence 
counting rate, with the counters and source straight line. The scattered 
coincidence rate will increase one counter moved around the 90° posi- 
tion which will lead errors the observed angular correlation. The photons 
scattered the backward direction the case the gamma rays have 
energies about 200 kev. this case possible discriminate effectively 
against the pulses from the scattered radiation and still allow most the 
pulses due the gamma rays (1.17 and 1.33 Mev.) counted. Accordingly, 
the discriminators were set eliminate pulses corresponding particles with 
energies less than 200 kev. The absence scattered coincidence counts was 
demonstrated placing the counters the 90° position and comparing the 
rates with and without lead shield between them. 


Readings were taken 10° intervals ranging from 90° 180°. The effi- 
ciencies the counters may vary slightly, owing variation the photo- 
multiplier high tension supply. Since the first approximation the coincidence 
rate proportional the product the two gamma counting rates, the slight 
variations efficiency may corrected for dividing the coincidence 
counting rate this product. For counters subtending very small solid 
angle the source, the coincidence rate function proportional 
the angular correlation function W(@) the relative probability 
emission the second photon angle with respect the direction 
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emission the first photon. However, the counters must subtend finite solid 
angle, necessary apply correction the coincidence rate order 
obtain W(@). The correction was calculated from the geometry the coun- 
suggested the uncorrected curve and consistent with the theory. The 
angle subtended the gamma counters was measured, using annihilation 
radiation from source. included angle 20° was selected order 
give sufficiently high counting rates without too great loss resolution. 
the case Co® the correction the measured correlation coefficient 
+0.147 was +0.006 giving 


The disintegration Sc“ follows scheme similar that the cascade 
gammas being only slightly lower energy. The techniques discussed for 
therefore, apply well The results the measurements Co® 
and are plotted Fig. The curve shown was calculated from 
theory (5) for two successive quadrupole transitions involving nuclear states 
with angular momenta: These results are good agreement with 
those obtained Brady and Deutsch (2). 


II. Beta-Gamma-Angular Correlation 


Since the beta particles from disintegrating nucleus have continuous 
energy spectrum, might expected that the corre- 
lation (if any) would function beta particle energy. The theory de- 
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Fic. gamma-gamma-angular correlation functions for and The solid 
curve calculated from the theory mentioned the text. 
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veloped Falkoff and Uhlenbeck (4) indicates that, general, this the 
case. Consequently, have distinguish between the differential angular 
correlation function and the integrated angular correlation function 
E). The function (6, represents the relative probability emission 
gamma photon direction with respect the direction emission the 
preceding beta particle, small range beta energy energy the 
beta detector sensitive all beta particles energy greater than then 


Eo 


where the number beta particles the energy interval (E, 
and the maximum energy the beta spectrum. 


The integrated angular correlation was measured using beta counter con- 
sisting layer anthracene flakes approximately mgm. per cm.? thick- 
ness RCA-5819 photomultiplier. The counter was practically 100% 
efficient beta counter but had low sensitivity gamma photons. The 
gamma counter was the same used previously the gamma-gamma cor- 
relation experiments. The source the form sodium antimonate was de- 
posited from solution zapon film, using insuline obtain uniform 
thickness, and was about 0.5 mgm. per thickness. was enclosed 
diameter vacuum chamber with the beta counter placed immediately 
outside aluminum window mgm. per cm.? thickness. 


First beta-gamma angular correlation experiment was performed with 
test for the presence scattering. The distribution was isotropic within 
the statistical accuracy (which was 1%) agreement with Novey (8). 


Next the integrated beta-gamma-angular correlation function for was 
measured. The decay scheme complex, involving five beta groups 
(9). However, the beta particles below Mev. are absorbed, the beta-gamma 
coincidence rate will involve mainly the 2.37 Mev. group with admixture 
coincidences due the 1.6 Mev. group; 270 mgm. per aluminum 
was used absorber front the beta counter, with which arrangement 
beta particles having energy greater than about 0.9 Mev. were counted. 
Readings were taken intervals 10° from 90° 180°. The gamma-gamma 
coincidence rate, which appeared background amounting 11% the 
beta-gamma coincidence rate, isotropic and has been subtracted. The cor- 
rection for the finite angular resolution 20° has been applied and the results 
plotted Fig. The vertical lines indicate the statistical deviation. The solid 
curve represents: W(@) 0.23 cos? which essential agreement with 
Beyster and Wiedenbeck (1) but somewhat different from the results ob- 
tained Ridgway (10). 


The differential beta-gamma angular correlation coefficient 


was measured the following method. anthracene crystal in. 
thick with RCA-5819 photomultiplier was employed beta counter 
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Fic. integrated beta-gamma-angular correlation function for for the beta energy 
larger than Mev. The solid curve represents the function: 1-0.23 
give pulses proportional the energy the beta particles. The work 
Hopkins (6) indicates that, experimental setup similar this one, 
the pulses are proportional beta energy least 3.2 Mev. The 
counter bias was set discriminate against pulses corresponding beta ener- 
gies less than about Mev. This bias also discriminates effectively against most 
the gamma pulses; however, gamma-gamma coincidences amounted 
this case 15-20% the beta-gamma coincidences. The beta pulses were 
transmitted through gating circuit channel kicksorter Chalk River 
design (7). The gating circuit was controlled the output the coincidence 
mixer, that was opened whenever coincidence occurred. The beta pulse 
which was delayed 0.25 usec. allow the gate open, was then transmitted 
the kicksorter. The kicksorter channels were evenly spaced over the 
upper part the beta spectrum from 1.0 Mev. 2.37 Mev. The accidental 
coincidence counts each kicksorter channel were obtained measuring the 
shape the beta pulse spectrum with the kicksorter and calculating from the 
total beta count which was obtained during each run. 


| 
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The photomultiplier high tension (1000 v.) was stabilized with battery 
reference and was maintained within 0.2 This was sufficient variation, 
however, cause slight change the beta pulse distribution, coincidence 
readings were taken alternately for hr. periods the 90° and 180° positions, 
the effects any drift present would minimized. The results the four 
individual runs agreed the worst cases within two standard deviations. 
(See our previous communication this connection (3).) The coincidence 
rates were divided the product the single channel rates compensate 
for any change counter efficiencies due changing the position the 
gamma counter. correction was also applied for decay the source. The 
gamma-gamma coincidence rate this case was not quite isotropic, being 
about 10% greater the position 90° than 180°. This effect was due 
scattering from one counter the other, the gamma counter necessarily being 
sensitive low energy gamma photons. The background gamma-gamma 
coincidence rate was subtracted give the true beta-gamma correlation co- 
efficient and the correction was applied for the angular resolution 20°. 


The energy scale terms beta pulse size was fixed determining the 
end point the beta pulse spectrum obtained the kicksorter and con- 
structing trial Kurie plot for each the four runs separately, the maximum 
beta energy being taken 2.37 Mev. According the results Langer, 
quoted (12), the beta spectrum (presumably the most energetic 
group it) the ‘‘a-type’’. Consequently, when constructing the Kurie 
plots, corresponding forbidden shape correction factor was introduced, al- 
though within the accuracy the kicksorter measurements this hardly 
necessary. The Kurie plots obtained, one which shown Fig. could 
resolved into two approximately straight lines, the lower energy beta group 
being agreement with the accepted value 1.6 Mev. The accuracy the 
end point (2.37 Mev.) the high energy beta group was about 0.05 Mev. 


The results the four separate series measurements, each involving four 
five days’ counting, have been averaged, and the curve for obtained 
plotted Fig. (The points corresponding the separate measurements are 
also indicated the figure.) check the consistency the measure- 
ments, the curve for and the beta spectrum determined the kick- 
sorter were used calculate (by numerical integration) the integrated angular 
correlation coefficient for 0.96 Mev. The value obtained, —0.26, 
satisfactory agreement with the value —0.23, measured directly, re- 
ported above. may also mentioned that the values reported 
Stevenson (11) for two values are consistent with the curve for 


shown Fig. 


view the decay scheme, lack accurate experimental data con- 
cerning the shape the beta spectra involved and the low accuracy the 
curve, unambiguous comparison this curve with theory seems hardly 
possible, this stage. are believe that the most energetic beta spec- 
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575 775 975 
Mev.- 1.56 1.97 2.37 


Fic. The Kurie plot the beta spectrum obtained with the kicksorter. The scale labelled 
represents the kicksorter bias arbitrary units. The energy scale Mev. proportional 
being adjusted Mev. the end point the the number counts 
correction factor. the total energy units the maximum energy the same units, 
and the momentum units mc. 


trum a-type (see remark above), than not unreasonable assume 
that the theoretical curve above 1.6 Mev. (end point the less energetic 
beta spectrum) would mainly determined the matrix element the 
Fermi theory. this assumption, and taking the angular momentum quantum 
number the ground state (even-even nucleus) zero, one still 
has variety possibilities regarding the angular momenta the other two 
involved. Fig. the two dotted curves represent the theoretical 
curves reported previously (3), which correspond two tentative assignments 
angular momenta: and course the curves may 
have significance (if any) only for energies above 1.6 Mev. All other assign- 
ments angular momenta consistent with the above assumptions seem 
definitely contradicted our experimental results; least, the gamma 
radiation involved multipolarity not higher than two. 
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Endpoint 
second beta 
spectrum 


COEFFICIENT 


1.0 1.2 1.4 1.6 2.0 2.2 
BETA Mev. 


Fic. beta-gamma angular correlation coefficient a(E) function the beta energy 
Experimental points corresponding different measurements are indicated symbols 
The average value a(E) with the standard deviation indicated the symbol 


the solid curve. dashed curves are calculated from the theory mentioned the text. 
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THE MAGNETIC DIPOLE OVER THE HORIZONTALLY 
STRATIFIED EARTH! 


Abstract 


The behavior small current-carrying wire loop over horizontally stratified 
earth investigated. The layers are considered have contrast conducti- 
vity and dielectric constant only. Both harmonic steady-state and step-func- 
tion current sources are considered. 


Introduction 


The earth’s crust example electrically dissipative medium since 
possesses finite conductivity. The conductivity exists virtue the result- 
ing total effect the electronic and the electrolytic current conduction the 
application electric field the medium (7, For many purposes 
this effective conductivity constant with respect frequency over wide 
limits. 

Depending the nature the soil rock the conductivity may vary from 
rather extreme high onemhoper rather extreme low value 
per meter. The conductivity may sometimes estimated from geological 
maps based the previous results tests for similar geologic formations. 
study these maps may indicate what type variation resistivity with depth 
expected. For example, there top layer sand with shales 
greater depth, the conductivity likely increase from very low value near 
the surface high value for greater depths. the other hand the surface 
layer clay while the underlying rock granite the opposite type variation 
can expected. The connection between earth conductivity and geological age 
and formation has been investigated Card (1) who examines the results 
large number measurements. 


The relative dielectric constant varies from low value two for dry rock 
clay high value for moist soil clay. The results measurements 
made radio frequencies for various rocks and soils are given Heiland (7, 
The displacement currents are usually very small for all significant 
frequencies compared the conduction currents. For this reason the conduc- 
tivity changes will far overshadow permittivity changes least for power and 
geophysical frequencies (i.e. frequencies less than one megacycle). 


The relative permiability rock formations and soil generally within 
unity for nearly all cases (7, 312). Except for static magnetic fields (i.e. 
the earth’s )the permeability always appears product with the conductivity 
(i.e. and hence again the conductivity variations greatly overweigh perme- 
ability changes (10). 


Manuscript received May 1961. 
Contribution from the Department Electrical Engineering, University Toronto, 
Toronto, Ont. 
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For many cases the earth can represented horizontally stratified 
medium with homogeneous and isotropic properties each layer. The analysis 
direct current flow such medium requires the solution Laplace’s equa- 
tion. The conditions the boundary and the current source sink are 
sufficient determine the potential everywhere. The classical solution this 
problem for two layer earth has been solved Maxwell (11) who used infinite 
series image sources with modified moments. Ollendorf (12) and Stefanescu and 
Schlumberger (19) have given solutions the form integrals involving Bessel 
functions. Formulas terms images have been given Hummel (8) for 
three layers. also gives approximate curves for this case. similar analysis 
given Ehrenburg and Watson (3) for any number layers. 


maximum information obtained concerning the subsurface conductivity 
variation when spectrum frequencies used for the source current. The 
behavior the observed field low frequencies characteristic potential 
theory. the other hand, the behavior.at high frequencies due partly 
propagation effects and characteristic wave theory. 


When source time varying current situated the vicinity layered 
earth the field components U(t) must satisfy the wave equation 
2 
the media, the boundary conditions the interfaces, and behave the proper 
manner near the current source and infinity. 


Sommerfeld (18) has derived general expressions for infinitesimal electric and 
magnetic dipoles energized harmonic current the surface semi-infinite 
flat homogeneous earth. The case the electric dipole has been extended the 
two layer case Riordan and Sunde (16) Foster (4) has developed mutual 
impedance formulas for finite wires the surface homogeneous earth 
double integration Sommerfeld’s horizontal electric dipole. Formulas for the 
field due transient current infinite wire the surface flat homo- 
geneous earth have been published Ollendorf (12) and Peterson (14). Riordan 
(15) has investigated the corresponding case for parallel finite wires. West (23), 
using numeral method, has evaluated the steady state impedance for certain 
special case two collinear grounded wires the flat homogeneous earth. 

this paper the formal solutions for harmonic time dependence will 
carried out for the vertical magnetic dipole over the two Interesting 
cases will then treated suitable specializations and approximations the 
above. The transient solutions will also given for several cases. 


Magnetic Dipole Two Layer Earth 


magnetic dipole small loop oriented with its axis the direction and 
located the origin the cylindrical coordinate system 2). 

The medium above the upper plane interface will characterized con- 
second units will used throughout. 
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The vector fields and everywhere can written terms single 
magnetic vector potential defined Schelkunoff (17, 128)* and are 
given 


The time factor exp will omitted. 


The primary field the same that due magnetic dipole strength 
IdA where the loop circulating current and the area the loop, 
Since the primary exciting current has only component then correspondingly 
the magnetic vector potential has only component. The primary vector 
potential then given Schelkunoff (17, 338), and has only component 
the direction. 


and 


Now the resulting magnetic vectors the three media must 
solution the wave equation: 


Solutions for these equations follow from the theory cylindrical harmonics. 
They are made solutions the following type: 


cos 
sin 


the Bessel functions the second type Y,(Ap) are not permissible. There 
obvious axial symmetry the problem that the plus minus sign the 
exponential must chosen reject the solutions that are not finite 
tends plus minus infinity. 


The general expressions for the magnetic vector potentials are then obtained 
integration over all values These vectors have only component the 
direction. 


The function occurring here related simply the often used magnetic Hertz potential 
are connected follows: 


2c 
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The characteristic functions (i.e. eigenfunctions) will now determined from 
the boundary conditions which require that the tangential electric and magnetic 
fields are continuous across the plane interfaces. using Equations (1) and 
(2) the tangential fields are given as: 


aF, 
The equations imposed the boundary conditions are then: 


These conditions hold for all values hence they can integrated with respect 
The function and their derivatives vanish tends infinity that 
the constant integration must zero. The conditions then reduce the 
following four equations; 


apply these conditions the primary excitation must expanded terms 
infinite integral involving cylindrical wave functions. This representation 
makes use the Sommerfeld integral (18) 


(15) 


Sommerfeld used this representation but did not prove for 
can checked however shown Appendix 


The four equations (11) and (12) then can solved for (A), (A), and 
The result this algebraic process is: 


| 
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the surface the earth (i.e. the magnetic vector potential then 
written: 


where the circulating current the loop and its area and 


R(A) = Alu, + Me wn 


The electric field also the surface the earth shown Fig. 
given by: 
‘i 


The voltage induced the element then equal E,ds. 


view the magnetic dipole and the receiving linear wire element the surface 
the earth. 
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the principle reciprocity given Schelkunoff (17, 478) the voltage 
shown Fig. given by: 


da x 


Fic. Plan view the electric current element and the receiving loop. 


Stokes’ theorem this induced voltage also given 


and 


The and terms are the components undetermined solenoidal field. 


The mutual impedance then between the current element and finite 
circuit between the points and along contour shown Fig. given 


Fic. The scheme for formulating mutual inductance wires the surface the earth. 
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where two dimensional vector whose components are and 
unit vector the direction the wire element the wire contour 


When and are the same point that the contour closed loop the 
latter integral zero since solenoidal field. The mutual impedance 
then given 


The above equation defines the mutual inductance between the two circuits. 


The mutual inductance between two coplanar circuits length and 
with angle between their directions the surface the earth then given 
by: 


The elements and are considered both one part closed loop. 


The mutual inductance between two coplanar circuits contour and 
the surface horizontally stratified earth then written 


The general formulation will now specialized various ways facilitate 
evaluation. 


Homogeneous Flat Earth 


The upper and lower layer are become identical that and 
constant that The for the magnetic vector poten- 


tial anywhere within the earth then reduces to: 


That is, all displacement currents are neglected. 


This can written the following form: 


| | 
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The integral involving has been already encountered and the integral involving 
Nhas been derived Foster (4) who carried out similar analysis for the electric 
dipole the homogeneous earth. Wolf (24) has also derived this latter integral. 


When the vector potential becomes 


The following relations are for the 
and (30) 


The required differentiations and then yield the relatively simple expres- 
sion for F,: 
Qn 
The electric field the surface the earth has only component and 
given 


(31) 


IdA 
The mutual impedance between wire element length and wire loop 


area the surface homogeneous earth then given by: 


TO, 
where the angle the direction the wire element and the line joining 
the wire element and the wire loop. 


Proceeding similar manner that used the case the electric dipole 
the infinite medium carried out the author previous paper (21), the step 
function mutual impedance (i.e. voltage response one circuit unit step 
function current the other) given by: 

where some small positive real constant. 


exp (34) 


Integrations are carried out and the results combined. The step function 
impedance given by: 


A(t) sin u(t), (35) 


and integrals are evaluated Appendix The character- 
istic function plotted Fig. 


WAIT: DIPOLE OVER HORIZONTALLY STRATIFIED EARTH 585 


(t) 


Fic. The step function mutual impedance between small linear wire element and small wire 
loop the surface homogeneous conducting earth. 


The vertical magnetic field the surface the earth can also found 


Now satisfies the wave equation this can written: 


The voltage e(w) induced another small loop area dA, distance from the 
source loop area then given by: 
Then the steady state mutual impedance between two horizontal loops area 
and dA, the surface the homogeneous earth given by: 
dA,dA 

The step function again given the contour integral, Equa- 
tion integrations are then carried out using the integrals listed 
Appendix 


and 
The characteristic function P(t) plotted Fig. 


| 
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Fic. The step function mutual impedance between two small wire loops the surface 
homogeneous conducting earth. 


Poorly Conducting Upper Layer 


many cases the conductivity the lower layer will much higher than the 
conductivity the upper layer There will then range 
frequencies where propagation effects the upper medium are negligible (i.e. 
This assumption will then also valid for the air 
regions above the lower conducting layer. effect this simplification Equation 
(8) rewritten with the following substitutions: 


Then the magnetic vector potential between the surface the earth and the 
surface the conducting layer given 


The field components anywhere the range can then found: 


(43) 


H,=0, 


possible represent terms tabulated functions valid for fairly 
wide range frequency and conductivity values the lower layer. The ex- 
pression (41) for the magnetic vector potential for this case can written 
the following way: 


586 
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ab (44) 


These three integrals (I), and (III) are then respectively given 


rao 


(46) 

where the second integral given Equations (29) and (31), 


2 2i 
(47) 


has been tabulated Foster (4) connection with mutual impedance finite 
grounded wires over homogeneous earth. 


the surface the earth then the magnetic vector potential given 
by: 
where the second term correction due the fact that magnetic dipole 
height above the conducting plane. This expression should compared 
Equation (31). From Equation (24) seen that mutual inductance 
between two coplanar circuits length and for this case 


Dipole over Thin Conducting Sheet Infinite Extent 

similar manner that used the problem the magnetic dipole over 
two layer earth the fields dipole over three layer earth can calculated 
cylindrical coordinate system used with the axis again perpendicular 
the 


= 
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The conductivity all regions outside the finite layer thickness depth 


and 


The term square brackets the latter integral can written follows: 


Now the case will considered where the thickness the conducting layer 
becomes very small. Also the conductivity becomes very large. The follow- 
ing expansions are then made: 


Now allowed approach zero and approach infinity such way 
that the product finite then the square bracket term becomes simply: 


da 


these limit processes considered finite parameter since the total 
integrand converges very rapidly for large values any case. 


The magnetic vector potential can then written: 


A+u 
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The second term Equation (52) represents the secondary field. con- 
venient transform the integral the following substitution: 


integration with respect can immediately carried out yield: 


This agreement with formula given Anton Graf (6) who used rather 
different approach. gives several expansion formulae for the secondary fields 
which are useful for computation for this steady state case. 


similar analysis has been carried out Stefanescu (20) for the case 
horizontal electric dipole over thin conducting plate. His method similar 
the method used here for the magnetic dipole. gives his fields the first 
order 1/g only and makes the usual assumption that the surrounding medium 
has negligible propagation characteristics. 


The transient fields will now calculated for the case when the loop ener- 
gized step function current. 


The steady state magnetic vector potential can written the following 
form: 


where 


The transient vector potential response f(t) for step function current 
then given by: 


The change the order integration here permissible that integrating 
first with respect there results: 


(57) 
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The integration with respect can now easily the unit 
impulse function defined Gardner and Barnes (5). 


The electric field component given by: 

(59) 

The step function impedance A(t) between loop area and 


current element length the surface the earth where there conduct- 
ing plate depth given 


where the angle between the direction and the direction the line 
joining the center the loop and the center the wire element. The 
characteristic function Z(t/hp, p/h) plotted Fig. and given 


A(t) 
(60) 


step function mutual impedance between wire element and wire loop above 
conducting sheet for various values the ratio p/h. 
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Here the first term characterized unit impulse function for the 
case when the medium infinite resistivity. Actually, owing the finite 
conductivity the medium surrounding the conducting sheet, the primary 
impulse will have finite duration. Its nature can expected similar 
the step function mutual impedance between the loop and the wire element 
the surface homogeneous earth finite conductivity given Equation 
(35). any case appears that the measuring time where 
the conductivity the surrounding earth rock material that the predomin- 
ant contribution the field due only the eddy current term involving 


Z(t/hp, p/h). 


then noted, observing the curves Fig. that for given conducting 
sheet depth the maximum transient response the point where the ratio 
p/hisunity. also seen that the decay time the field considerably 
greater for larger values the ratio p/h. 


Appendix 


Since exp solution the wave equation, finite infinity, and 
has axial symmetry, then must the form 


the Fourier Bessel theorem 


Therefore 


The free space case when purely imaginary special case 
integral which given Watson (22, 416). 
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Appendix 
well known contour integral given Churchill (2) his tables Laplace 
transforms: 


where and where erf the error function argument Each 
side this equation can differentiated three times, with respect yield 
the following three integrals: 


271 c-i@ 


where erf’ and are the first and second derivatives, respectively, the 
error function. The error function and its derivatives are defined and tabulated 
Jahnke and Emde (9). 
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THE QUANTIZATION THE CLASSICAL THEORY 
SPINNING 


SHANMUGADHASAN? 


Abstract 


The classical theory particles, possessing charge and dipole moment, and 
moving electromagnetic field, considered the assumptions that there 
constraint connection between the rotational variables and the velocity 
the particle, and that the two invariant squares the dipole moment six—vector 
are constants the motion. different schemes are obtained according 
the two invariant scalar products the dipole moment and total spin angular 
momentum are are not constants the motion. The 
Corben theory fits into the former scheme. The classical schemes are put into 
canonical form using for each particle the relativistic connection between 
the momenta and the rest-mass, modified include the effect the kinetic and 
potential energies due spin and dipole moment, the Hamilton—Jacobi 
equation and the usual Poisson brackets for the translational and total spin 
variables. The Wentzel field and the \-limiting process are used mainly deal- 
ing with the field. variational principle for the equations 
given with the field treated according the limiting process Dirac the 
relativistic cutoff method Feynman. The quantization completed 
using the analogy rules. The changes required when the interacting field 
vector meson field are 


Introduction 


The quantum theory dynamical system obtained first putting the 
classical equations motion the system into the canonical (Hamiltonian) 
form, and then quantizing the theory using the standard rules connecting 
the canonical form classical quantum dynamics. Since quan- 
tum mechanics the dynamical variables the system are linear operators satis- 
fying noncommutative algebra while classical mechanics they satisfy 
commutative algebra, this procedure for quantization may lead ambiguities 
the order which the factors must appear the quantum theory. spite 
this difficulty the method can used obtain definite results many 
problems. 


The relativistic quantum mechanics spinless electrons obtained the 
above method, the elementary Lorentz classical theory the electron yielding the 
Gordon-Klein wave equation the quantum theory. But the spinning electron 
occurring nature described correctly the Dirac linear wave equation for 
particles with spin half quantum, and not the Gordon-Klein equation. 
view the success the Dirac equation little interest has been shown ex- 
tending the Gordon-Klein equation include the effect spin and dipole 
moment. The first attempts setting relativistic classical theory spin- 
ning particles possessing dipole moment were made Frenkel and 
Thomas. But there has been systematic formulation the corresponding 
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quantum theoretical equations. therefore worthwhile from the point 
view theoretical completeness investigate such equations. 

Following Dirac his theory radiating point electrons, Bhabha and Corben 
(2) have set general classical theory spinning particles proceeding from 
the assumptions that the dipole moment proportional the spin the par- 
ticle, and that the invariant square the spin tensor constant during the 
motion. Without making use these assumptions the present author (10) has 
arrived general form for the rotational and translational equations motion 
when there constraint relation between the rotational and translational 
variables. This scheme not sufficient determine the motion the particles 
without further assumptions regards the rotational motion the particles. 
Such assumptions can made referring situations the three dimensional 
dynamics tops and gyroscopes. pure gyroscope the magnitude the 
total spin angular momentum three-vector remains constant. The relativistic 
analogue this has been dealt with (10). The present paper deals with 
rotational motions more complicated and possessing more rotational degrees 
freedom than pure gyroscopic motion. Then the total spin angular momentum 
vector represented the sum the so-called spin and extra vector 
angular momentum. This extra vector can chosen many ways. But two 
simple cases are obtained choosing this vector that (i) its magnitude 
constant and the angle between the two vectors variable, (ii) its magnitude 
variable and the angle between the two vectors constant (particularly with 
the two vectors orthogonal). The relativistic analogues these two cases are 
considered below. Since the second case corresponds symmetric top, 
leads the Bhabha-Corben equations. 

The classical theory can put into canonical form the equations motion 
follow from variational principle (3, then possible introduce canon- 
ical coordinates and momenta, and hence obtain the Hamilton-Jacobi equations 
and the Poisson brackets (P.B.’s hereafter) the dynamical variables. How- 
ever action principle not necessary, provided that the suitable Hamilton- 
Jacobi equations and consistent scheme P.B.’s can postulated from general 
considerations (4, 5), the cases considered below. Yet interest 
have the action principle for least one case, and the Bhabha-Corben scheme 
chosen for this purpose. Quasi-coordinates appear the covariant form 
the action principle. 

The present-day forms electrodynamics are not entirely satisfactory, be- 
cause difficulties arise dealing with the interaction the particle with its own 
field. For the sake definiteness the discussion the present paper confined 
mainly the form electrodynamics using the Wentzel field and the 
process. There are features the theory considered below that are unaffected 
these questions, and they could taken over into satisfactory form field 
dynamics. Modifications the treatment the field have been suggested 
recently several authors. The changes the action principle necessary 
incorporate the limiting procedure Dirac (6) the relativistic cutoff procedure 
Feynman (7, are briefly indicated. 
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The extension the theory when the interacting wave field (Proca) vector 
meson field quite straightforward most respects, and has been dealt with 
the P.B. relations for the field variables given there are not correct, 
because they lead inconsistencies when dealing with the extra condition in- 
volving these variables. This difficulty resolved here postulating suitable 
for the field. 


the quantum theory the Schrédinger equations for the particles are natural 
extensions the Gordon-Klein equation. The physical interpretation there- 
fore must based the Gordon-Klein probability rule which not positive 
definite. the wave equations the field intensity tensor appears coupled with 
the spin. This not the linear wave equation for electron unless 
modified, Pauli has done, give the particle abnormal magnetic moment. 
the present theory the same type wave equation applicable particles 
any spin and any value the dipole moment. 


remains seen whether the theory applicable any the elementary 
particles known experimentally. The quantum theory analogous pure 
gyroscopic rotation the classical theory not suitable for electron, has 
been seen from the computation the Compton effect total cross 
section for Compton scattering increases indefinitely with the energy the inci- 
dent photons, which also true the cross section based the Pauli modified 
linear wave equation. Since the heavy particles like the nucleons have ex- 
plicit spin interaction with the meson field, the theory may interest this 
connection. The calculation scattering cross sections, which can used 
check the new quantum schemes with the facts nature, yet done. 


The Classical Equations Motion 


Relativistic notation used most the following work. departure from 
the usual convention tensor analysis the one (adopted Dirac this type 
work) denoting contravariant quantities subscripts, and covariant 
quantities superscripts. The velocity light taken unity, that 
the metric tensor (Greek indices take values has all components 
used whenever convenient. point space-time then denoted 
The square written The notation for the scalar product 
and (a, b). convenient some stages the work introduce 
three-vector notation. three-vector denoted sans-serif bold type. Thus 
are used for the scalar product and vector product respectively two three- 
vectors 


With tensors the second rank, and are denoted and 
(XY) respectively. The antisymmetric tensor formed transvection from 
two antisymmetric tensors X,,, Y,,, namely denoted 
The antisymmetric tensor dual (some writers call supplement- 
ary) the antisymmetric tensor defined terms X,, 
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where the permutation symbol antisymmetric each pair indices and 
such that the above definition, follows that X,,. 


Further, X,,, Y,, are two antisymmetric tensors, then 


Any given antisymmetric tensor gives rise two self-dual tensors 
(the factor being introduced for convenience) 


since 
self-dual tensor has three independent components only, since 


The three independent components, the (23), (31), (12), the (01), (02), (03) 
components form three-vector. through this step that three-vectors 
where Y,, antisymmetric. Then 


equation the type where X,, and Y,, are above, can 
written terms self-dual tensors the form 


Consider typical particle having coordinates z(s) and velocity (with 
Here the proper-time, and the dot denotes differentiation with 
respect tos. Let the charge the particle and let its dipole moment (both 
electric and magnetic dipoles) described the antisymmetric tensor Z,,, 
which may considered function The total spin angular momentum 
the particle described the antisymmetric tensor field quanti- 
ties are obtained from the four-potentials according F,, 
where denotes The actual field any point expressed terms 
the ingoing and retarded fields, the outgoing and advanced fields, the 


When the system consists several particles, the label affixed the quan- 
tities associated with the 7th particle. 


The starting point the present work the general form, given (10), for 
the translational and rotational equations motion particles the above 
type, when constraint connection between the rotational variables and the 
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velocity the particle assumed. The translational equation motion has 
the form 


where scalar representing the total mechanical mass the particle, and 
the rotational equation motion, the form 


The effective field f,, the above equations motion single particle the 
mean the ingoing and outgoing fields. When there are several particles, the 
equations motion the ith particle have exactly the same form above, 
except that the effective field now 


(3) 
The invariant equations that can deduced from (1) and (2) are 
dm/ds (4) 
(NZ) =0, (Nf) (5) 


Two definite schemes equations motion corresponding the two cases 
mentioned the introduction are obtained below. 


Case Suppose that Z,, CS,,, where constant, and that 


The variation R,, with consistent with the constancy and (RR*), 
given 


where X,, some antisymmetric tensor yet determined (Ref. (1)). Then 

Since (SS) Equation (7) leads 


now assumed that independent R,,, and that R,, independent 
Then, for (8) true for R,,, the condition 


must hold every point the world-line the particle. 


The six homogeneous linear equations equivalent (9) can written the 
matrix form 
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The antisymmetric matrix singular. Since its rank must even, 

most four. The determinant the lower right, bounded the third row 

and third column, and the last row and last column, has the value 

which not zero general. Thus the rank the matrix four. The 

solution the system equations therefore contains two arbitrary parameters. 

X,, S,, evidently solution the system. Since the antisymmetric tensor 

the left side (9) vanishes, follows that the tensor dual this also must 

vanish. Hence [S*. X],, which gives X,, The general 
solution the system (9) 


where are two arbitrary scalar parameters. 


view (10), Equation (7) becomes 


result, and (RR) (RR*) (RS) (R*S) (4) leads 


If, grounds simplicity, the arbitrary parameters are chosen constants, 
constant. The constants and are the nature moments inertia. 
The second and third terms the expression for represent the kinetic energy 
due rotation the particle, which naturally appears relativistic theory 
addition the constant rest-mass the particle. The translational 
motion given 


while the rotational motion given the two equations 
Case II. Assume that N,, and Z,, satisfy the relations 
(NZ) constant, constant, (15) 
Note that and (NZ*) are not constant Case 
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where Y,, some antisymmetric tensor, assumed independent Z,,, and 
which yet determined. Since (NZ) from (5) and (15), (17) gives 
=0. For this hold everywhere the world-line for 
arbitrary Z,,, Y],, must true all points the world-line. Hence 
Y,, DN,, where and are arbitrary scalars. For simplicity, 
and are chosen constants the work that follows. 


With the above choice Y,,, Equation (4) becomes 


constant. and are again analogous moments inertia, 
and the rotational kinetic energy represented the second and third terms 
the expression for appears addition the rest-mass The motion 
the particle thus described 


for the translation, and 


for the rotation. closer examination these equations motion, with the 
help self-dual tensors, shows that can expressed terms Z,, and 
and their duals. Then the scheme equations reduces the same form 
that Bhabha and Corben (2). 


Introduce now the three-vectors arising from the self-dual tensors. The 
corresponding sans-serif bold lower case letter denotes the three-vector from 
given six-vector; for example, corresponds self-dual tensor with the 
like U,, used here, though there particular reason for 
this preference. Then Equation (20) becomes 


Hence 
that is, 
Also, 


These three-vector equations can immediately translated into tensor equa- 


tions. Equation (21) corresponds 
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Hence has the form 


Similarly, from (22), 


that and (NN*) have the form 


terms and (ZZ*), (ZN), and which are related such 
way that the equations motion are consistent. Bhabha and Corben have 
used exactly the form (23) for the total spin angular momentum and the ob- 
tained using (24) the total mechanical mass the particle. 


The Bhabha-Corben equations can considered from the following point 
view. Put Equation (2) Z,, CS,,, where constant, and S,, the 
spin angular momentum about the axis spin. The total spin angular momen- 
tum symmetrical top ordinary dynamics (for instance Ref. (9)) 
K[s where the so-called spin along the axis (with constant), 
and the moment inertia perpendicular the axis. The relativistic 
analogue this S,, where S,, the spin angular momentum 
along the axis spin (with constant), and constant. this value 
chosen for N,, Equation (2), the rotational equation motion becomes 


with the invariant equation K(SS). Hence 
translational equation becomes 
given terms N,, and S,, 


the expression for one could add constant multiples the dual tensors 
S*, and with the subsequent modifications give the general form 
the Bhabha-Corben equations. But simpler consider the form (25, 
26) the Bhabha-Corben equations. 


The classical equations motion considered far apply one particle only. 
The generalization system several particles quite straightforward. The 
typical ith particle has the equations motion exactly the same form 
those for single particle except that the label added all the constants 


Zz 
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and tensors describing the ith particle, and that the effective field given 
Equation (3). 
The effective field all the above equations motion can expressed 


terms new field, called the Wentzel field. The Wentzel potential, 
particles possessing charge and dipole moment Z,,; has been defined Ref. 
(10) 


where A(x) the Jordan-Pauli relativistic defined 


A(x) 2x,6(x’). 


(27) 


then obtained from the potential 

where small time-like four-vector made tend zero ultimately, pro- 
vided that 


(29) 


The potential defined (27) and (28) satisfies the equations 


(30) 


Canonical Form the Classical Theory 


The general (canonical) equation motion dynamical variable 
system F], where denotes the P.B. and and where 
the Hamilton-Jacobi equation the system. change from the in- 
dependent variable new independent variable would necessitate the use 
Fds/dr the new Hamilton-Jacobi equation. order apply this 
result particular problem, one must know the Hamilton-Jacobi equations 
the system, and the P.B. relations connecting the various dynamical variables 
describing the system. According Dirac (5), permissible postulate the 
P.B.’s without the explicit introduction canonical coordinates and momenta, 
provided that the postulated relations satisfy the following properties 
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Case Consider the equations motion (12)—(14). These equations reduce 
those spinless electron when R,, and S,, are zero for all The required 
Hamilton-Jacobi equation must such that reduces that spinless elec- 
tron, which known the Gordon-Klein type, when R,, and S,, are put 
constant plays the role total effective mass the particle (12), 
conformity with relativistic ideas that the kinetic energy spin and the poten- 
tial energy due the dipole moment the field appear additions the rest- 
mass the spinless particle. Hence the required Hamilton-Jacobi equation 
taken that obtained from the Gordon-Klein type replacing the constant 
rest-mass the effective mass the particle, namely, 


where the momentum four-vector, and A(x) the field, Equation (28). 


The choice the P.B.’s not difficult. The usual relations for the co- 
ordinates and momenta, and those for the spin angular momentum variables are 
taken over. One notes also that the translational and rotational variables 
the particle refer different degrees freedom. Thus for the particle variables 


[z,, spin variables] [p,, spin variables] (36) 

R,, has P.B. relations similar those From the P.B.’s for S,,, one gets 


Similar relations hold for and for 


(37) 


Since the particle variables and the field variables belong different degrees 
freedom, A,(x) has vanishing with z,, S.s, The P.B. two field 


where D(x, x’, defined (33). Further discussion the field equations 


deferred the end this section, since they have the same form both the 
cases considered here. 


The correct connection between velocity and momentum given 


that 
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(40) 


Equations (39) and (40) together give the correct translational equation 
motion 


The equations motion for S,, and R,, are 


which describe correctly the rotational motion the particle. Also, S’, (SS*), 
R’, have vanishing P.B.’s with all the dynamical variables the system, 
and therefore with which should hold for them constants the motion. 


Case The Hamilton-Jacobi equation suitable for dealing with Equations 
(18)—(20) 


with the field A(x) given (28). and have the usual P.B. relations (35) 
among themselves and have vanishing P.B.’s with and The P.B. 
relations for are taken the usual ones for spin angular momentum. 
taken play the role generalized coordinates relation the angular 
momentum Thus the for the rotational variables are 


and the corresponding relations involving the dual tensors. The field variables 


again have the same P.B.’s Case Then, since have vanishing 
P.B.’s with (ZZ*), (NZ), and and since 


the equations the motion check out correctly, and (NZ), (NZ*) 
are constants the motion. 


The equivalence the equations under consideration the Bhabha-Corben 
equations has already been shown. Since the for Z,, can obtained from 
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(42) using (20), the above discussion covers the canonical formulation the 
Bhabha-Corben equations. However, separate discussion the canonical 
formulation Equations (25) and (26) given below. 
Introduce the three-vectors and connected the relation 2K[s Xs]. 
Hence 
The P.B.’s satisfied and are 
[Say Sy] 0, So] Sey = Ney (44) 


where the suffixes take values that are cyclic permutation 
Therefore, from (43) and (44), has the P.B.’s 

where the usual Kronecker symbol. terms the original tensors, the 
P.B.’s for the spin variables are those corresponding (42), and 

(46) 


where 


and A,(x) have the same P.B.’s before. The above scheme P.B.’s 
and the Hamilton-Jacobi equation 


lead the canonical formulation Equations (25) and (26). 


The generalization the case several particles can carried out the 
same lines the work Dirac, Ref. Any dynamical variable 
general function all the proper-times the particles. The equations 
the particle. For the consistency these equations, the condition 
F;] must satisfied. The correct equations motion are obtained, 
assumed that (i) the P.B. two variables referring different particles 
zero, (ii) the P.B. the field variable A,(x) and any particle variable zero, 
(iii) the P.B.’s the variables referring any particle are the same form 
for particle, and (iv) the Hamilton-Jacobi equation for each 
particle the same form for single particle. The dynamical variables 
(here they are tensor quantities) the particle depend only and are 
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assumptions. The equation just gives the equations motion 
the particle, since they are exactly the same form for single particle. 
Also easily verified that the consistency condition [F;, F;] satisfied 
when the inequalities (29) hold. 


The equation motion for the field variables dA,(x)/ds; [A,(x), F;]. 
This agrees with Equation (32), provided that [A,(x), A,(x’)] g,,D(x, x’, 
result which has been postulated earlier. Hence 


[A.(X), fur (x’)] (- LarPy + £..0,) D(x, d), \ (48) 
[fas(X), fur (X’)] = (- + 2.9.9, + — £5,0,0,)D (x, x’, 
Equations (30) and (31), which are not consequences the equations 
motion, have imposed extra conditions. They are consistent with the 
equations motion, since, Ref. (10), 


[OA, F;] = 0, [R(x), F;] = 0. 
Equation (31) plays the role the Hamilton-Jacobi equation for the field. 


There one these equations for each and they are consistent with each 
other since [R(x), R(x’)] 


Formulation the Action Principle 


The action principle for the translational motion spinless electrons well 
known (for instance (3)). the case particles with spin and dipole moment, 
the action principle must give both the translational and rotational equations 
motion. Before proceeding the solution the general problem, simpler 
consider first the situation when the particle has only rotational motion with- 
out translation. should then possible combine this action function with 
the action function for the translational motion, the form which known from 
the action function for spinless electrons, give the complete action principle. 


terms the three-vectors, the Bhabha-Corben rotational equation 
motion (25) takes the form 


Before discussing (49), consider first the action principle for the equation 
motion 


Since assumed that there translational motion, the dot Equations 
(49) and (50) denotes differentiation with respect the time. 

Let Oxyz fixed orthogonal triad ordinary three-space. moving triad 
having Eulerian angles obtained from position Oxyz rotating success- 
ively through about Oz, through about the new position Oy, and through 
about the new position Oz. The formal extension complex Eulerian angles 
complex three-space (analogous the one above) used below the 
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further discussion. the fixed complex three-space the velocity the 
representation Since constant, has the form 


when the total angular velocity given 


being unit vectors along the directions about which the rotations 
are respectively performed. Hence (s, Note that 
arbitrary (52). 

Now treat the quasi-velocity vector corresponding quasi-coordinate 
vector the true coordinates the system being that 


The and not commute when operating quasi-coordinates 


The transitivity equation for the quasi-coordinate vector is, with the help 
(53) and (54), 


Consider the action integral 


previous certain time and form its variation allowing vary. Since 
from (55), 

and from (51), 


one gets, from (56), 


where denotes the variation denotes the total variation the final 
integral sign (57) gives the correct equation motion (50). 


The terms under the integral sign cancel when the equation motion used. 
The first term the remaining ones the right hand side (57) perfect 
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differential, and can therefore discarded. the standard procedure 
for introducing canonical coordinates and momenta, the momentum con- 
jugate the canonical coordinate and the energy the system. 
Hence [¢, s,] result that can used evaluate the P.B.’s the variables 


which are just the P.B.’s postulated Ref. (10) for the discussion Equation 


(50). 
Consider now Equation (49). Take the action integral 


and perform the usual variation before. Then 


The coefficient under the integral sign (59) gives the correct equation 
motion. 


With the help (51) and (53), (59) becomes 


The first term the right discarded, since perfect differential. The 
other terms give the canonical coordinates Q,, and their respective momenta 


Since 
the P.B.’s the various dynamical variables can evaluated from (60) and 
(61). Thus 


{P, cos Q,} sin cosec Q,] 
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{P, Q,} cos cosec Q,, sin Q,] 

cos cos cosec cos sin Q,} 


and similar relations, which together lead the P.B.’s (44) and (45) considered 
earlier. 


The real action function for the rotational motion known immediately 
combining the three-vector complex form and its complex conjugate. The action 
function for both the translational and rotational equations motion ob- 
tained generalizing the action function for the translational motion spinless 
electron include the action function for the rotational motion. However, 
must remembered that, dealing with the action function for rotational 
motion, the equations have been considered ordinary Euclidean space whose 
underlying metric the form dx; With such metric there 
difference between covariant and contravariant quantities. The translational 
and rotational equations have been worked out Lorentz space, where the 
scalar product two four-vectors has negative sign the terms involving the 
space components. But the metric the three-space used 
then the value the rotational action function this space can taken 
over immediately into Lorentz space. The change from the ordinary Euclidean 
space the new pseudo-Euclidean space means replacing every typical scalar 
product m), say, m). This will then automatically show the differ- 
ence (in sign) between the contravariant and covariant components any 
vector. Thus the action function for the rotational motion has taken 
over The form the P.B.’s, however, remains the same before. 


The field equations can also obtained from the action principle. The usual 
infinities associated with point model the particle have been avoided the 
above work using the Wentzel field and the process. The action 
principle for this case can given the same manner Dirac 
ever, more general form limiting process due Dirac (6) can used 
the action principle (62). 


parametric problems the calculus variations, the integrals varied 
must independent the parametric representations used. This condition 
satisfied the integrands are functions which are positively homogeneous 
degree one the derivatives with respect the parameter. the discussion 
the action principle here, the action integrals for the particles involve the 
proper-time parameter. The homogeneity condition can satisfied for the 
action integral introducing suitable powers factors terms where 


| 

| 

( 
7 

= 
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this step called for, done below. Then the action integral proper 
parametric form, and change any parameter other than the proper-time 
possible. 


The action integral for the equations motion several particles and the 
field 


where dx,dx,dx,dx,, and the denotes summation over the several par- 
ticles (the label omitted for convenience). (62) P,, the quasi-velocity 
corresponding the quasi-coordinate Q,,, that dQ,, P,,ds. Further, P,, 
connected with the equation 


other properties for are assumed further work. The action function for 
the field, namely, the last integral (62), has the Fermi form, thus avoiding the 
zero value for the momentum conjugate A,. The details 
handling the action principle with the above type field equations are found 
Ref. (6). 


The relativistic cutoff procedure Feynman (7, can taken over into 
the present theory considering the action integral 


where f,, the field tensor derived from the potential four-vector 


behaves like for large dimensions but differs for small, and 
normalized that 


The properties with special examples, and their consequences have been 
discussed Feynman. 


Quantization 


passing the quantum theory all the classical dynamical variables are 
replaced linear operators, with real classical variables corresponding self- 
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adjoint operators. The operators operate vector |>, called ket-vector, 
many-dimensional space, which vector specifies the states the dynamical 
system quantum mechanics. The coordinates any representation 
that can set give the Schrédinger wave functions. The linear operators 
are general noncommutative, and satisfy commutation relations correspond- 
ing the classical P.B. relations according the formula Ba) 
This noncommutative character the operators gives rise ambiguity, 
whenever there product two factors, which have nonvanishing P.B. 
the classical theory and therefore noncommuting the quantum theory. 
simple examples, considerations simplicity are invoked choose the order 
the factors leading consistent theory, the example considered below. 
The Hamilton-Jacobi equations and supplementary conditions give operator 
equations called the Schrédinger wave equations. All the other equa- 
tions become relations between the operators the quantum theory. When 
there are several wave equations the type they are consistent 
provided that the F;] for all 


When the dynamical system consists several particles and the wave field, 
the above rules lead wave equations the form 


(64) 


for each particle, where has one other the forms (34), (41), (47). The 
supplementary condition (31) gives the wave equation 


R(x) (65) 


for the field. The extra condition goes over unchanged into the quan- 
tum theory relation between the corresponding operators. The wave 
equations are all consistent, since has been verified that 


the quantum theory corresponding the Bhabha-Corben classical equa- 
tions, (25) and (26), linear operators corresponding their classical counter- 
parts are the dynamical variables satisfying commutation 
relations corresponding (35), (46), (38), and (48). more convenient, 
however, use and (or and their commutation relations correspond- 
ing (44) and (45) set the matrix representation the spin variables. 
quantum theory, the order the factors appearing shown. From the com- 
wherever occurs the classical theory. If, instead the 
process, the limiting process corresponding the action principle (62) used, 
then would necessary write the wave equation for the field terms the 
Fourier components the field variables exactly the same way Ref. (6). 


The relativistic cutoff procedure for quantum electrodynamics corresponding 
the action principle (63) can given the same lines Ref. (8). 


AS 


SHANMUGADHASAN: QUANTUM THEORY SPINNING PARTICLES 


Extension the Theory Vector Meson Field 


The extension the preceding type theory the case when the interacting 
field vector meson field satisfying the equation has been 
considered Ref. (10). The classical equations motion, the Wentzel field, 
and the Hamilton-Jacobi equations have the same form the electromagnetic 
case except that the vector meson field variables U,, replace the variables 
fy», respectively, and the generalized Jordan-Pauli function A*(x), defined 


where the Bessel function order one, replaces the A(x). 


the classical theory, the electromagnetic field equations (with can 
obtained from the meson field equations (with letting 
This not true, however, the quantized field equations. known that 
the P.B. relations for the electromagnetic field cannot obtained from those 
for the meson field letting Further, the P.B.’s for the meson field 
can chosen that they are consistent with the condition analogous the 
Lorentz condition, while this not the case with the electromagnetic field, 
where the Lorentz condition has introduced supplementary condition. 
Ref. (10), the P.B.’s postulated for the field variables and the subsequent 
treatment the extra condition 


correct treatment these questions follows. 

Split the potential U,(x) the form 


where 
(68) 


Since 

Wentzel potentials, 


But 


Equations (69) and (70) agree provided that the P.B. relations for *U, are 


taken 


Hence 
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(72) 
and therefore 


Also, from (71), with *G,, 


[*Gas(X), *G,,(x’)] (— — 


Equation (72) shows that the extra condition (66) can taken over unchanged 
into the quantum theory. the Hamilton-Jacobi equations, and therefore 
*G,,. The for *U, and *G,, are given (71) and (73). 


The Stueckelberg variation the above formulation split *U, the 
form 
(74) 
equivalent 
(75) 
The P.B. relations for and are taken 


Equation (75) not consistent with (76). But since 


Equation (75) treated supplementary condition (analogous the electro- 
magnetic case). 
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